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SUMMARY 

The work statemeiit items of this contract have been addressed in 

this^Bi^ report in the following manner; 

A) A thorough bibliography and summary of significant results has 
been presented in Appendix A of thislSlBiBi report. 

•B) A discussion of the more significant problems associated with 
digital computer generation of holograms has been presented 
in Chapters I and II. 

C) A criteria for producing optimum digital holograms has been 
presented in Chapter II. This criteria revolves around amplitude 
resolution and spatial frequency limitations induced by the 
computer and plotter process. 

D) Test results have been compiled and examples of the products 
desired have been created digitally and compound with optically 
produced products. Chapter IV presents these results. 
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CHAPTER I 


INTRODUCTION 

In classical holography, a recording is made of the interference 
pattern of an object wave with a reference wave. This system requires 
the use of a laser, an optical bench, and some recording medium. The 
binary computer generated hologram proposed by Lohmann and Paris^^* 
presents several advantages over the classical -hologram: (1) a physi- 

cal object need not be used, (2) no laser is required and (3) the 
requirement for expensive optical equipment is eliminated. In realiz- 
ing computer generated holograms, however, physical limitations of the 
computer must be considered. 

Storage capacity of a digital computer limits the. number of terms 
in the frequency domain representation of the binary hologram. The 
device which plots the binary mask for making the hologram is con- 
strained to operate in discrete steps. Thus, the dynamic range of the 
Fourier transform numbers which are represented in the binary hologram 
is limited by the digital computer and its associated input-output 
devices. The problem of producing a binary hologram becomes one of 
satisfying the physical limitations of the digital computer while 
maintaining the best possible image reconstruction. 

One facet of binary holography that has not been exploited widely 
is its use in cryptography; that is, the encoding of information in 
binary holograms.® A cursory study of encoding/ decoding binary 


*Superscripts apply to entries in the bibliography of Appendix A. 
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holograms has revealed numerous difficulties in retrieving the encoded 
information. The alignment of the spatial filters and the related 
optical equipment is particularly laborious. The problems encountered 
with the optical system naturally lead to the idea of simulating the 
optical system with a digital computer as a means of studying infor- 
mation degradation and as a way of optimizing image reconstruction. 

The fact that lasers, optical benches, and recording media are not 
required makes computer simulated holography particularly appealing. 

The computer simulation process to be used in this study is 
analagous to a coherent optical processing system^ capable of 
realizing operations of the form' 


I(x,y) 


K 



g(C,n)h(x 


2 


C»y - n) d? dn 


( 1 ) 


where I(x,y) is an intensity distribution, K is a complex constant, 
g(x,y) is the object function, and h(x,y) is the system impulse 
response. 

A configuration for realizing the intensity distribution I(x,y) 
is shown in Figure 1 . Source S represents a coherent light source and 
lens serves as a collimator. The input to be processed is inserted 
as a space varying amplitude transmittance g(x,y) in plane P2. (In 
the computer simulation, the binary object mask appears at this 
point). Lens L2 Fourier transforms g(x,y) producing an amplitude 
distribution 


KjG(x2/lf,y^Af) 


( 2 ) 



'3 . 



Figure 1 . Configuration for, Optical Processing System, 

at P^, the back focal plane of is a complex constant, 

G = F[g], and X is the wavelength of the source. 

The amplitude and phase of the frequency domain distribution at 
may be manipulated by the- insertion of a frequency plane filter . 
with amplitude transmittance 


t(x2,y2> = KgHCx^/Xf .y^/Xf) 


( 3 ) 


where H = F[h]. Thus the amplitude distribution to the right of P^ is 
simply the frequency domain product GH. The lens Lg Fourier trans- 
forms GH to yield an intensity distribution at P3 


I(x3,y3) 


K 


CO 



g(C>n)h(-x3 


2 


C,-y3 - n)dc dn 


( 4 ) 


where the minus signs preceding X3 and y3 are a consequence of 


FlFlg(x,y)JJ 


g(-x,~y) . 


( 5 ) 
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In particular, if the transmittance of the frequency domain filter at 
Is t(x 2 ,y 2 ) = 1» then GH = G and lens Lg simply reconstructs the 
original image ftmction from with the space coordinates reversed. 

•This sign problem is resolved easily by defining the coordinate 
system and y^ as shown in Figure 1. 

A particular application of the optical processing system of 
Figure 1 is in the area of character recognition. A specific 
spatial signal s(x,y) is inserted into the system at Pj^. A matched 
filter is then inserted at P^.- By definition, a linear space invariant 
filter is said to be matched to a signal s(x,y) if its impulse 
response h(x,y). is given by h(x,y) = s*(-x,-y). For an input signal 
g(x,y) applied to a matched filter h(x,y) , the output v(x,y) is 


^(x,y) =11 

— OO 

g(C,n)h(x - 5 ,y - n)d? dn 

(6) 

CO 

v(x,y) =11 

g(s>n)s*(? - X, n - y) d 5 dn 

(7) 


oo 


which is just the cross correlation of g and s. Fourier transforma- 
tion of the impulse response of the matched filter requires that 



F[h(x,y)l 

= F[s*(-x,-y)] 

(8) 

or 

H(f ,f ) 

X y 

= S*(f ,f ) . 

X y 

(9) 


Note that denotes the complex conjugate. Thus the frequency plane 
mask at P^ must have amplitude transmittance proportional to S*. The 
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field distribution transmitted by the mask is then proportional to 



s(c,p)s*(? - X, n - y)dc dn 


SS* 


( 10 ) 


which is entirely real. This implies that the frequency plane mask 
cancels all curvature of the wave front yielding a plane wave. This 
plane wave front is then brought to a bright focus at by the lens 
L^. Thus the presence of a signal may be detected by measuring the 
intensity of light as a function of input signal s. Optical 

interpretation of the matched filtering operation is shown in Figure 2 
and the output of the processor for a cross correlation application is 
shown in Figure 3. and are the spatial width of h and g respec- 
tively. 


Computer Generated Holograms 

Generally , computer , generated holograms are capable of functioning 
just as optically formed holograms. The following is a discussion of 
some of the basic principles of computer generated holograms. 

The Whit taker-shannon sampling theorem states that it is possible 
to compute a limited number of samples of a continuous function and 
with these samples to reconstruct the continuous function exactly. 

The preceding theorem is subject to the constraints that (1) the con- 
tinuous function is band- limited and (2) the continuous function is 
sampled at least twice in any increment as large as the spatial period 
of the highest spatial frequency in the function. Thus, it is possible 
to obtain a Fourier transfom representation of an appropriately 



BlXiy) transmittance 

S*(x,y) 


Figure 2. Optical Interpretation of the Matched Filtering 
Operation. 





Figure 3. Location of Various terms of the Processor 
Output with a Matched Filter. 
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sampled object function. In realizing a computer generated hologram, 
the computer must calculate a great many samples of the two-dimensional 
Fourier transform of the subject. In making binary holograms by means 
of a computer, the transmittance is considered to be zero or one. The 
transmittance of the binairy mask representing the subject is also one 
or zero. 

In order to understand how the computer can create a hologram 
which is comparable to an optically produced hologram, it is necessary 
to consider the mathematics of sampling and of the discrete Fourier 
transformation. Mathematically, sampling may be accomplished by means 
of the comb function,' defined as 

oo 

comb(x) = y 6(x - n) (11) 

where 6(x) is the unit impulse function and n is an integer. See 
Figure 4. 


i 



coml 

i 

>(x) 

. i 



[0 1 2 


Figure 4. The Comb Function 
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Furthermore,' the Fourier transform of the coirib function, 

F[comb(x)] = comb(f) . ■ (12) 

The sampling operation may be represented as a multiplication of a 
continuous function g(x) by the function 

“ 1 

^ 6(x - nAx) =“ — comb(x/Ax) . (13) 

AX 

n=*-^ 

where Ax is the sampling interval. Thus, the sampled function 
g^(x) may be expressed as ' 

g (x) = g(x).-^ comb(x/Ax) 

S AX 

CO 

= g(x) I 6(x - nAx) 
n=-~ 

oo 

= I g(nAx), 6(x - nAx) (14) 

n=^ 

Each sample of g(x) is a delta fimction whose strength is given by the 
value of g(x) at the position of the delta function as shown in 
Figure 5. , • 


AMPLITUDE 

4 



Mx— H 


Figure 5. The Continuous Function g(x) and the Sampled 
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Computer generated holograms require computation of the Fourier 

transform G (x) of the sampled spatial function g (x) . If the sampling 
s s 

has been carried out in accordance with the sampling theorem, the 

Fourier transform G(f) may be obtained from G (f) . The original 

s 

spatial function is simply the desired image reconstruction and may be 
accomplished optically or mathematically, the latter being the 
approach used in the proposed simulation process.- 

In transforming 8g(x), the product g(x) and ^ comb(x/Ax) in the 
spatial domain becomes the convolution of their Fourier transforms in 
the frequency domain. 

G (f) - G(f)*comb(Axf) 

s 

CO 

= . G(f)* ”5! 

n=-» 

From the definition of convolutionj 


G (f) 
s 

G^(£) 


oo 

— 5 ; 

Ax ^ • 


G(u)6(f “ n/Ax - u)du 




Ax 


' CO 

i ^ G(f - n/Ax) , 


n=- 


(16) 

(17) 


as shown in Figure 6. 


ATX Z 


G^(£) 


G(f)- 


x:^ ATX 


‘•max 


max 


Ax 


H 


Figure 6. Fourier Transform G (f) of the Sampled 
Spatial Function. 
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If the original space function g(s) is bandlimited so that. G(f) 
has non- zero values only in the interval 


£max. £ ^ ^max 

2 - “2 


( 18 ) 


overlap of the shifted transform is prevented as long as ^ max 

is satisfied. Thus, ix < 1/f specifies the conditions of the 

— max 

Whittaker-Shannon sampling theorem. 

To recover G(f) from G (f ) , .the sampled Fourier transform is 

simply multiplied by a rectangular window function Ax ’ 

\ 

G(f) = Gg(f)Ax rect(f/f ^^^ ) . '(19) 


In order to recover the original spatial function g(x), it is 
necessary to take the inverse Fourier transform of the product 
G^(f)Ax ^ect(f/f^^). This is equivalent to convolving in the spatial 
domain the inverse Fourier transforms of these two functions* The 
inverse Fourier transform of the window function is given by 


F ^[Ax rect(f/f )] 
max 


Axf 




max -Iff X 
max 


( 20 ) 


Thus, 


g(x) = F“l[6(f)J = gg(x)* Axf^ 


sin(7ff x) 
max 

irf X 
max 


^ g(nAx)6(x - nAx) * Axf ^^.^ 


sin (irf x) 


max 


n=-“ 


irf , X 
max 


00 sin irf (x “ nAx) 

r, . ^ max 

^ f fe - nto) 




( 21 ) 
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The preceding discussion may be extended to two dimensional spatial 
functions g(x,y) by replacing the comb function with a two dimen- 
sional array of delta functions, coinb(x) comb(y). 

In practice, the Fast Fourier transform algorithm is used for 
computing the Fourier transform of the object function. In the 
production of computer generated holograms, the Fourier transformed 
function representation is plotted and photographed to create a 
binary Fourier transform hologram. It is worth noting that the 
binary hologram is a special form of spatial filter. 



CHAPTER II 


■ DEFINITION OF THE PROBLEJ 

The primary sources of problems in generating synthetic holo- 
grams' are the physical limitations of the digital computer and its 
peripheral equipment. The main objectives of this research are to 
(1) study the effects of computer induced degradations on binary 
holograms, (2) .develop a technique for optimizing reconstruction of 
images by computer generated holograms, and (3) verify the optimization 
process by observing image reconstructions, both optically and by 
computer simulation. 

Some knowledge of the Lohmann-type hologram is necessary in 
order that the problems of generating holograms by means of the com- 
puter be fully appreciated. Lohmann holograms are binary in nature 
and have a non-negative, real amplitude transmittance.^^ The 
Fraunhofer' effect is achieved by means of the so called "detour 
phas.e." The synthetic hologram produced by this- procedure is in 
effect a diffraction grating which will reproduce an image when 
illuminated by a coherent, monochromatic light source. Image recon- 
struction from. a Fraunhofer hologram is illustrated in Figure 7. The 
Fraunhofer diffraction pattern in the image plane is described by the 
function 

00 f X f y 

f , _j2 ( -^ + -2- ) 

U(x,y) = J J F(f^,fy)e df^df^ ; (22) 




Binary 

Hologram 


Reconstructed 

Image 


Figure 7. Optical Setup for Reconstruction of Binary Hologram. 
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X represents the wavelength of the incident light and f is the focal 
length of the lens. The problem becomes that of synthesizing the 
complex function 


F(f , f ) = 

X y 


f X f y 

U(x,y)e d^dy 


(23) 


which may be recognized as the inverse of the two-dimensional Fourier 
transform of equation (22). Lohmann and Paris have presented a 
method of synthesizing binary synthetic Fraunhofer holograms capable 
of realizing the image U(x,y). Of course, there are some constraints 
placed on the image, limiting both its size and niimber of frequency 
terms . 

Limitations on image size require a finite extent Ax = Ay in 
area. The minimum resolvable region is limited to 6x = 6y = 1/Ax. 

The number of resolvable points in the image is then given by 

= St = (f >' • ' 

The niimber is known as the space-bandwidth product and is also 
equal to the number of cells in the binary hologram. Each cell of the 
Lohmann hologram acts as a miniature diffraction grating; its geometry 
is shown in Figure 8. 

Since both the object and the reproduced image are finite in 
extent, the two-dimensional Fourier transform of. equation (22) may be 
replaced by a two-dimensional discrete Fourier transform .of terms. 
The object to be described is divided into N x N regular grids with 
spacing 6v = 1/Ax, Each of these space samples will have a complex 
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frequency term F(ra6v,n6v). In turn, this complex number will be 
recorded as dictated by Lobmann’s algorithm to produce a synthetic 
binary hologram. 

The complex frequency domain number representation is recorded 
in the Lohmann cell of Figure 8. The magnitude of the transmittance 
is determined by the size of the cell aperture, . In practice, 

this is the only magnitude parameter allowed to vary, C6v remaining 
fixed at approximately 1/2 the cell size for maximum image intensity. 
The phase of the complex frequency term determines the position of the 
aperture, . The cell can shift from -tt to fr, modulo 2ir, with the 

center of the cell representing zero phase. 

At Mississippi State University programs have been developed for 
generating Lohmann-type holograms using a Univac 1106 computer and a - 
Gould electrostatic plotter. Cell realization is based on the ability 
of the plotter to plot 80 dots per inch in two dimensions. The field 
chosen for the synthetic hologram is 10 inches by 10 inches, thus the 
the plotter provides an array of 800 by 800 dots for each hologram 
produced. The cell realized by the Gould plotter is shown in Figure 9 
The physical extent of each Lohmann cell is limited by the number 
of frequency temns in the hologram. In this study, the sample space 
is 64 by 64 and the corresponding maxlmimi number of frequency terms 
is the same. Thus, a cell size may vary from 800/64 = 12 dots to a 
maximum size of 800 dots for a hologram containing only one frequency 
term. Again, it should be noted that the cells are square by design. 
Also observe that in recording the 64 frequency terms in 800 dots the 
number of' dots per cell is 12 due to the plotter quantization. 
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Figure 8. Lohmaim Cell Geometry. 





Figure 9. Gould Plot of Lohmann Cell. 


POOB QUALIIY 
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The smallest change that can be made in the amplitude represen- 
tation in the case of 64 frequency terms is 1/12 = 8.33% quantization 
in amplitude. For a hologram recording with only 40 frequency terms, 
the cell size will be 800/40 = 20 dots resulting in a quantization 
error of only 5%. Thus, the smaller the number of frequency terms, 
the larger the cell and the less the quantization error. 

Phase quantization may be regarded in a similar manner. In the 
example of 64 frequency terms and 12 dots per cell, the phase quanti- 
zation is 2Tr/12 = tt/6 or 30® phase quantization steps. Reduction to 
40 frequency terms and the corresponding 20 dots per cell yields phase 
quantization of 2tt/,20 = 18®. Again, larger cells reduce the quanti- 
zation error. 

Generating the hologram becomes a problem in terms of the 
dynamic range of the frequency terms. The higher frequency terms are 
usually small in amplitude, thus much of their information content is 
lost in quantization steps . One technique for reducing dynamic range 
is to band limit the number of frequency terms. In this method care 
must be exercised not to sacrifice information contained in the higher 
frequency terms for the sake of Improved resolution. 

A second method of reducing dynamic range of the hologram is to 
amplitude limit the low frequency terms of the Fourier transform 
numbers. Since a binary mask resembles an array of rectangle functions, 
the Fourier transform of each point is basically a sine function. The 
zeroiil order term is usually large in amplitude with the higher fre- 
quency terms decreasing rapidly. The low frequency, large amplitude 
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terms may be clipped to reduce the dynamic range of the hologram 
recording without appreciable loss of information. 

In order to obtain the best possible image reconstruction 
within the limits of the digital computer, some image quality 
measurement must be defined. For this investigation, an integral of 
signal amplitude squared over the reconstruction aperture has been 
chosen. This technique relates directly to intensity error and 
is easily adaptable to digital computer analysis. 

This error measurement entails a comparison of the pure recon- 
struction image to the image aperture which has been clipped and 
band-limited. The algorithm for determining the numerical error 
introduced by clipping and bandlimiting first Fourier transforms the 
object aperture under investigation; the resulting frequency domain 

function g(f ,f ) is again Fourier transformed to produce I(x,y), 

X y 

the space domain image with no degradation. The frequency domain 


expression g(f ,f ) is also operated on by amplitude clipping and 
X y 

bandlimiting to reduce its dynamic range. The corrupted array 

g (f f ) is then returned to the spatial domain as I (x,y). The two 
c X 5 y c 

resulting arrays are then term by term compared, squared, and summed , 
over the image aperture to yield the error figure. A block diagram 
of the error algorithm is shown in Figure 10. 

Since the error criterion is based on image intensity, only the 
amplitude of the image apertures need to be considered. To avoid the 
possibility of scale differences in the two reconstructed images being 
compared, a scaling factor has been used in the error analysis 
algorithm. The discrete error function including the scale factor s 
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Figure 10. Block Diagram of Error Analysis Algorithm. 


may be written 

I I (|I„J - (25) 

m n 


where and are the pure image and the degraded image 

respectively over the discrete image aperture which is m by n. The 
problem of finding a minimum of the error function in terms of the 
parameter s is solved by expanding the squared quantity under the 
double summation, differentiating with respect to s, then setting the 
resulting expression to zero to solve for s. The resulting value of 
s which minimizes the error function is 


s 


n 

m n 


I • I 
^ mn cran' 


1 1 1 

m n 


cmn 


(26) 
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Finally, to obtain a signal to noise figure for the corrupted 
image, the error figure is divided into the pure image power 

nuj" ( 27 ) 

m n 

The signal to noise figure for a given set of clipping and band- 
limiting parameters is given by 

II II 1^ 

SN = - j , (28) 

y y (I - si 

" ^ mn cmn 

m n 

where the image is a function of the clipping and bandlimiting 
parameters chosen to satisfy computer constraints. 

One of the major efforts of this research is a computer program 
for determining an optimum set of parameters to be used in producing 
a synthetic hologram. The optimizing routine employs the signal to 
noise criterion~described above. The function to be optimized is a 
non-linear function of two variables with an additional constraint 
that the bandlimiting parameter be an integer number. Optimizing a 
function of this complexity is a formidable problem, and according to 
Dr. R. R. Hocking,* authority on non-linear programming techniques, 
"there is no guarantee that an absolute maximum or minimum of such a 
function can be located." A sequential search technique with a post- 
optimum search was eventually employed with success. 

4 

Simulated hologram reconstruction was implemented because of the 
difficulty of achieving good optical reconstructions with limited 

*Prlvate communication with Dr. R. R. Hocking, Professor of 
Mathematics and Statistics, Mississippi ■ State, MS, 39762. 




optical facilities. Simulated reconstructions allow images to be 
observed and signal to noise optimization routines to be studied 
without the usual photography, laser, and optical paraphernalia. 

The simulation program functions as a two lens optical system 
as shown in Figure 11. A flow chart for the simulation program is 
given in Figure 12. This routine may be used for various spatial 
filtering experiments as well as simple image reconstruction. For 
example, matched filtering and autocorrelation may be simulated by 
this algorithm. 

One of the major efforts of this research has been the recovery 
of the holograms from the tape which is used to drive the plotting 
equipment. The Information recorded on the tape has been clipped and 
bandlimited as previously described. In addition, it has suffered 
quantization effects due to the discreteness of the plotter. Thus, 
the recovery of this Fourier domain array makes available a simulated 
image with all of its degradation effects. The flow chart for this 
algorithm is in Figure 13. 
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Figure 11. Optical System with Two Lenses. 
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Figure 12. Computer Simulation of Two Lens Optical System. 
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Figure 13. Flow Chart of Reconstruction Program 







CHAPTER III 


A SYNOPSIS OF DIGITAL COMPUTER PROGRAMS 

The Univac 1106 Computer at the Mississippi State University 
Computing Center was used for the development of all holograms » 
spatial filters, and' simulation plots in- this study. Thus,, it is 
necessary to describe the programs and to explain their usage. The 
majority of the programs are written in Fortran V, the exceptions- 
being those programs directly related to the Gould plotter which are 
written in Univac 1106 assembler language. A listing of all programs 
and subroutines used in this study is given in Appendix B. 

There are some routines which are employed in virtually all 
computer runs; these .programs are on file and are called as sub- 
routines as required by a main program. The commonly used subroutines 
will be discussed- first. 

•The spatial dimension of the object description utilised in this 
study was chosen to be a 64 by 64 square array. This necessitates a 
complex dimensioning of 64 by 64 storage arrays in the computer. Com- 
plex specification, is required because the Fourier transformation of 
the binary transmittance masks produces complex. numbers. This study 
is restricted to two-dimensional arrays although it could be extended 
to include three-dimensional cases. Because the primary concern in 
this study is information storage and processing, the two-dimensional 
analysis is sufficient. 
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Subroutine FFT2(A,M) 

FFT2 is a two dimensional Fast Fourier transformation based on 
the Cooley-Tookey algorithm. The input matrix A is a M by M array of 
complex numbers in rectangular form. M may be any integer power of 
two up to a maximum value of 64. The ordering of the elements of A 
requires that A be operated on by subroutine ROTATE prior to insertion 
into FFT2. The matrix A is returned as the Fourier transformation of 
the input matrix. 

Subroutine ROTATE (A, B.M) 

ROTATE takes an M by M complex array A as read in row-column 
order of- a spatial coordinate system and rotates it into matrix B 
which is in the order required by FFT2. Henceforth, a reference to a 
rotated matrix implies that a matrix has been operated on by ROTATE 
and is in the order expected by FFT2. The input matrix A remains 
unchanged . 

Subroutine FOLAR(A.B,M) 

This subroutine takes a M by M complex array A which is in rec- 
tangular form (a + jb) and converts it .to the M by M matrix B in polar 
form. The magnitude of B is stored as the real component of B while 
the angle appears as the imaginary part of B, Matrix A is unchanged. 

Subroutine BUILD (A. B .NO .M) 

This program functions as a bandpass filter and limits the output 
matrix B to NT terms, where NT is less than or equal to 64. The output 
is matrix B, NT by NT, complex, and un-rotated. NT is simply the 
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desired bandpass and is a function of the input parameter NQ. NQ 
is related to the number of quanta in the Lohmann cell created by the 
plotter. NT is found from an array QT in subroutine BLOCK' and common 
QUAN to be described later. The input matrix A remains unchanged. 

The output matrix B is used by subroutine BITS to produce the binary 
hologram. 

Subroutine BITS(A.CLIP.NQ,H12,B) 

BITS takes the bandlimited matrix A from BUILD, and with a 
specified clipping level CLIP, forms the 800 by 800 bit array used to 
produce the binary hologram. The input matrix A must be in polar 
form. The array of terms is NT by NT, where NT is determined by 

■ ’ NT = QT(1,NQ) . (29) 

The cell size is found from the same data block and the parameter NQ 
by 

NC = QT(2,NQ) (30) 

The resulting phase quantization is 2tt/NC radians and the aperture 
width is set for optimum intensity as shown in Figure 8 to be NC/2 or 
C = 1/2. Note that the subroutine argument N12 is used to specify a 
positive or a negative type of plot; N12 = 1 plots a black background 
with transparent apertures while N12 = 2 causes the background to be 
white with opaque apertures. The output matrix B is simply a one- 
dimensional 4600 term array which is a portion of the bits used in 
making the hologram. 
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Subroutine INITm2,NC) 

This is an assembler' language routine which is called by BITS to 
initialize subroutine BS to a specified cell size NC and aperture 
width N2. . 

• Subroutine BS(B,IWL.IMAG,ISH) 

BS is another assembler language program used by BITS. Each time 
BITS calls BS one of the NT by NT bit cells of the hologram is gene- 
rated by bit manipulation. The argument list is specified by BITS 
and is of no concern to the user. 

Subroutine BITS2(A,CLIP.NQ.N12) 

This subroutine is used in forming the plot of the reconstructed 
image in the simulation program. The input matrix A should be in 
polar form. 

Subroutine CRUPT (B . C . CLIP .NQ .N) 

CRUPT takes the frequency domain matrix B, which is in polar 
rotated form and corrupts it as prescribed by parameters CLIP and NQ. 
The argimient N is 64 as set by the- maximum number of terms in the 
complete spatial array. The returned matrix C is NT by NT terms as 
determined by the argument NQ. C is returned in rectangular, rotated 
form. The matrix B remains unchanged. 

Subroutine OPT(A,B.F,NQ,PCT,X) 

This subroutine finds optimum clipping and bandllmitlng values 
of the object aperture to be made into a hologram. The matrix A is 
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64 by 64 and has been Fourier transformed twice. It is in polar 
form. The matrix B is in polar form, has been Fourier transformed 
once, and is 64 by 64. Matrix F is returned as an array of signal to 
noise figure while X is the optimum signal to noise figure. NQ and 
PCT are returned as optimum values for achieving the best signal to 
noise ratio in the s^thetic hologram. See Figure 14. 

Subroutine PIC(A,CLIP,NQ,N) 

This program plots a nomograph of the amplitude limited array 
which, was used to produce the hologram. The input matrix A is in 
polar, un-rotated form and is N by N in size. CLIP and NQ are used 
as in previous subroutines. 

Subroutine TAPE(L.$) 

This subroutine reads the bits from Fortran File 9 that are used 
in making the actual plot of the hologram. 

Function FMAX(A,N) 

The maximum value of a matrix A is found by this function sub- 
routine. The matrix A is N by N in size and in polar form. 

Function ERR(A.B.ASQ,CLIP.NQ,N) 

The error analysis previously described is accomplished by this 
program. Matrix A is the uncorrupted N by N reference image I in 
polar form. Matrix B is N by N, polar, and uncorrupted at its time 
of entry into ERR. CLIP and NQ apply as before to clip and bandlimit 
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Figure 14. Flow Chart for Optimization Program OPT. 
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Figure 14. (Continued) 
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Figure 14* (Continued) 









matrix B. ASQ is the reference image power that is used in the 
signal to noise figure calculation. 


Block Data BLOCK 

BLOCK initializes the named common QDAN which contains the QT 
array. It is from QT that the bandwidth NT and the cell size NC are 
found as given by equations (29) and (30) respectively. 

Main Program HOLPLO 

This program takes the 64 by 64 array of binary ntjmbers used to 
describe an object and creates a synthetic hologram. A flow chart 
for this program is given in Figure 15. HOLPLO makes use of several 
of the previously listed subroutines. 

Main Program SIMUL 

SIMUL is the Fortran program which simulates the two lens optical 
system. Its flow chart is given in Figure 12. The plotted output of 
SIMUL is the image reconstruction of the synthetic hologram created by 
HOLPLO. 


Main Program OPTl 

OPTl takes the object description which is to be made into a hologram 
and determine optimum clipping and bandwidth values for it. The outputs 
•NQ and PCT from stibroutine OPT may be inserted directly into HOLPLO for 
production of a near optimal hologram. ‘ 




Figure 15. Flow Chart of the Synthesis of a Computer 
Generated Hologram 
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Main Program 0PT2 

This program functions in the same manner as OPTl except that it 
finds the optimum value of an encoded hologram. 


Main Program SN 

Program SN varies the clipping levels and the number of frequency- 
terms in the Fourier transformed array to provide an output array of 
signal to noise figures as functions of variables NQ and PCT . Main 
program OPTl supercedes SN and eliminates the need for unnecessary 
and lengthy calculations. 

Main Program SECON 

BECON is used in conjunction with HOLPLO to recover the hologram 
from the tape which drives the plotter. In ordinary usage, RECON 
follows HOLPLO sequentially and provides a simulated image reconstruc- 
tion which includes clipping, bandlimiting, and phase and amplitude 
quantization induced by the plotter. See Figure 13 for a flow chart 
of RECON. 


Main Program SNPLOT 

SNPLOT is a routine for calculating the signal to noise ratio of 
the hologram after it has been degraded by the plotter. It is run 


sequentially after HOLPLO. 
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Haln Program SPAFIL 

SPAFIL is used for simulating an optical spatial filtering pro- 
cessor. The inputs are descriptions of objects which are compared 
with digitally generated spatial filters. The output is a plot of 
the simulated autocorrelation. This program may be used for other 
spatial filtering applications. 

Main Program CODE 

CODE is a program which modifies HOLPLO for creating an encoded 
computer generated hologram. Main program 0PT2 is employed to find 
optimum values for NQ and PCT. A SIMUL type program is used for 
recovering the encoded information and plotting the reconstruction. 



CHAPTER IV 

PRESENTATION OF RESULTS 

The synthetic holograms produced in this research effort are made 
by (1) electrostatically plotting the Lohmann-cells generated by Sub- 
routine BITS and (2) photographically reducing this 10 inch by 10 inch 
array to form the hologram. The physical size of the hologram is made , 
to be consistent with the optical bench used in the reconstruction 
process. Collimators, irises, and lenses are used as required. 

Main program HOLPLO is used to make the plot from which the holo- 
gram is formed. The photograph was made by Illuminating the plot with 
two 100 watt flood lamps.' The hologram was recorded on ortho-type 
film, exposure time one minute at / 5.6. The camera used was a 
Practina with standard lens. The exposure time and /-stop values were 
determined by .trial and error. Since the end product is a binary holo- 
gram, film non-linearities do not present a problem. 

For this report, the plotter outputs were chosen to have white 
backgrounds with opague apertures to minimize noise introduced by the 
plotting equipment. Figures 16, 17, and 18 show Gould plots for three 
different clipping and bandlimiting values. Figure 16 is the plot for 
near-optimum parameter values as determined by OPTl'. Figure 17 shows 
the plot set for the optimum number of frequency terms but with the 
amplitude clipped to 0.1 of the value of its maximum amplitude term. 

The plot of Figure 18. is for the case of no amplitude clipping but 
the number of frequency terms has been reduced to 18. 
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The object used in. making the previously described masks was 
the block letters AB . Optical reconstruction of the holograms pro- 
duced for the three cases was accomplished by use of the optical, 
arrangement shown in Figure 7. The reconstructed images were recorded 
on film located in the image plane. Figure 19 is the optical recon- 
struction with near-optimum parameter values. Figure 20 is the 
reconstruction with severe clipping of amplitude terms, while Figure 
,21 shows the bandlimited version of the reconstruction. 

Optimum parameter values for synthesizing holograms are found 
by use of main program SN . This program allows both the bandwidth 
and amplitude clipping values to be varied in order that the signal 
to noise ratio of the resulting hologram may be observed. Table 1 
shows signal to noise values of the hologram of AB as clipping and 
bandlimiting are -varied. The clipping level is set by truncating any 
magnitude term greater than some preset percent of the maximum ampli- 
tude term of the transform array. The level is set by 

PCX = N/9. , N = 1, 2, 9. (31) 

The number of frequency terms NT is determined as in Subroutine BUILD 
by parameter NQ, where 

NQ = 3*M , M = 1, 2, •••, 13. (32) 

The optimum clipping and bandlimiting parameters are determined by 
finding the maximum signal to noise figure from the chart of Table 1. 
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Figure 19. Optical Reconstruction of Near-Optimum 
Hologram. 
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Optical Reconstruction of Hologram with 
Severe Clipping. 
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Figure 21. 


Optical Reconstruction of Bandlimlted 
Hologram. 



Figure 22. Simulated Reconstruction of AB without 
Degradation. 



TABLE 1. Signal to Noise Values for the Hologram AB as 

Clipping and Bandlimitlng Parameters are Varied. 


PCT 

1/9 

2/9 

3/9 

4/9 

5/9 

6/9 

7/9 

■ 8/9 

9/9 

NQ 

3 

1.23133 

1.26116 

1.30343 

1.34114 

1 

1.37070 

1.39255 

1.40629 

1.41311 

1.41527 

6 

1. -37075 

• 1.42290 

1.47815 

* 1 

1.52640 

1.56493 

1.59325 

. 1.61160 

1.62049 

1.62085 

9 

1.53503 

1.63277 

1.70643 

1.77311 

1.82809 

1.87095 

1.90109 

1.91725 

i. 92186 

12 

1.73836 

1.88801 

1.98742 

2.07706 

2.15273 

2.20923 

2.25094 

2.27382 • 

2.27968 

15 

1.96549 

2.18778 

2.31967 

• 2.44189 • 

2.54602 

2.63231 

2.68817 

2.72577 

2.73226 

18 

2,21838 

2.51712 

2.68914 

2.85310 

2.99298 

3.10394 . 

'3.18864 

3.22791 

3.23327 

21 

2.49294 

2.90470 

3.13744 

3.35534 

3.54824 

, 3.71513 

3.80620 

3.86246 

3.85195 

24 

2.80401 

3.34233 

3.64422 

3.93939 

4.20581 

4.40654 

4.53874 ■ 

4.57331 

4.55456’ 

27 

3.16853 

3.88491 

4.28995 

4.67748 

5.03094 

5.27615 

. 5.38977 

5.41485 

5.31762 

30 

3.74277 

4.82634 

5.43430 

6.00485 

6.45327 

6.73444 

6.83529 

6.70898 

6.38879 

33 

4.27517 

5.73963 ■ 

6.48660 

7.17552 

7.54540 

7.72156 

7.55086 

7.13220 

6.68360 

36 

5.01664 

7.03466 

7.79711 

8.18133 

8.28672 

7.97956 

7.48325 

6.86054 

6.29099 

39 

5.42138 

7.26897 

7.27034 

7.11410 

6.68921 

6.25740 

5.82334 

5.30017 

4.83029 
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For the object employed in this study, the maximum signal to 
noise figure is 8.28672, as indicated in the table, and occurs 
when the clipping level is 5/9 of the maximum frequency domain 
number, and with NQ equal to 36. A value of 36 for NQ corresponds 
to 52 frequency terms, as determined by block data BLOCK and 
common QUAN. 

A more efficient method of finding the optimum signal to 
noise value for a S3aithetlc hologram and its associated clipping • 
and bandlimiting parameters makes use of the program OPTl and sub- 
routine OPT. OPTl simply reads in the binary mask from which the 
hologram is to be constructed. OPT is the non-linear optimization 
routine which performs the search for a maximum signal to noise 
value' as the two parameters PCT and NQ are varied. It should be 
noted that the signal to noise figure is a non-linear function of 
the parameters, and additionally that the parameters are constrained 
to vary in discrete steps. A computer listing of the signal to 
noise figures and the corresponding parameters as the optimization 
routine searches for a maximum is given in Table 2. The starting 
point of the search routine is arbitrary; however, based on some 
knowledge of the function and the desired output, the solution is 
found faster by starting with a high NQ value. The printed outputs 
are signal to noise values fo\and during the search, with the last 
line listing the maximum signal to noise figure along with the 
optimum parameter values . 



TABLE 2. Computer Printout. of Optimization 
Program OPT. 


NQ=39 

CL= 3 

S/N= 7,27034 

Na=36 

CL= 3 

S/Ns 

: 7,79711 

I=li 

U= 2 

S/N= 

5.73963 

1=11 

0= 3 

S/N= 

6.46660 

1=11 

0= 4 

S/N= 

7.17552 

1=12 

0= 2 

S/N= 

7.03466 

1=12 

J= 4 

S/N= 

8.18133 

1=11 

j= 3 

S/N= 

6.48660 

1=11 

j= 4 

S/N= 

7.17552 

1=11 

5 

S/N= 

7.54540 

1=12 

J= 3 

5/N= 

7.79711 

1=12 

J= 5 

S/N= 

6*28672 

1=11 

J= 4 

S/N= 

7.17552 

1=11 

J= 5 

S/N= 

7.54540 

1=11 

J= 6 

S/NS 

7.72156 

1=12 

J= 4 

S/N= 

6,18133 

11 

J= 6 

S/N= 

7.97956 

1=13 

J= 4 . 

S/N= 

7.11410 

1=13 

js 5 

S/N= 

6.66921 

1=13 

U= 6 

S/N= 

6.25740 


THE OPTIMUM BanDWIOTH AND CLIPPING LEVEL VaLUES ARE; 


Ha 

CLIP 

F(12» 5) 

36 

.5556 

8,28672 
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For comparison purposes, the three cases described above were 
reconstructed utilizing main program SIMUL, The first example simply 
illustrates the ability of the simulation routine to reconstruct an 
image without degradation. The uncorrupted reconstruction of AB is 
shown in Figure 22. The irregularities in the letters are due to 
quantization in the input binary mask. Figures 23, 24, and 25 are 
simulated reconstructions of optimum, severe clipping, and band- 
limiting cases, respectively. 

Main program EEGON is a simulation routine which recovers the 
hologram exactly as it was recorded by the plotter. Bandlimiting and 
clipping effects on hologram quality have been observed and studied 
in terms of a signal to noise ratio figure. RECON includes clipping 
and bandlimiting effects as well as phase and amplitude quantization 
introduced by the plotting equipment. Figure 26 shows a simulated 
reconstruction of AB after it has suffered plotter degradation. In 
addition, main program SNPLOT analyzes the signal to noise ratio of 
the degraded image; this will be discussed in the following chapter. 

In order to test the simulated optical processing system, the 

152 

famous Abbe-Porter wire mesh, experiment was performed. Figure 27 shows 
the binary mask used to describe the mesh and Figure 28 is a horizontal 
slit filter. The results are as recorded by the simulator in Figure 
29, If the slit is rotated by 90®, the results are as shown in 
Figure 30. 

Further experiments with matched filters were performed utilizing 
the simulated optical processor. The character T as shown in 
Figure 31 was used to produce a matched filter as described in 
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Chapter I. The simulated autocorrelation of an input character T 
with the matched filter T is plotted in Figure 32. The optical 
autocorrelation of T with the matched filter T is photographed in 
Figure 33. Observe that in the optical reconstruction there is a 
reimaging of the input character and the convolution of the two 
characters which do not occur in the simulated field of reconstruc- 
tion. Figure 34 is the result of the input character X which is 
not matched to the filter T . As in the previous example. Figure 35 
is a photograph of the optical autocorrelation of the unmatched 
filter and input. To further illustrate that the simulator functions 
as an optical processor. Figure 36 presents results of changing the 
position of the character T in the input aperture. 

The final experiment with matched filtering was with character 
recognition. Once again, the character T of Figure 31 was used as 
the filter. The mask OXT of Figure 37 was applied as the input. The 
simulated autocorrelation of the input characters OXT with the filter 
T is plotted in Figure 38. The darker plot on the right side of the 
figure aligns with the input symbol T and represents the auto- 
correlation of this input with the filter T . The center spot is the 
cross correlation of the input X with the filter T , and the left- 
most spot indicates the cross correlation of the input 0 with the 
filter T . An optical autocorrelation of the same character recog- 
nition system is photographed in Figure 39. 

The last effort in the study deals with encoding and decoding 
holograms. The characters AB are used as the message symbols. 

Object description AB is Fourier transformed by HOLPLO in the usual 
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manner. After the transformation, a reference symbol, which was 
chosen to be 0 for the study, is also Fourier transformed. A holo- 
gram is then produced of F[AB] divided by F[0]. The result is an 
encoded hologram F[AB]/F[0] which is plotted and photographed as 
previously described. Theoretically, the original information 
symbols are recoverable by optically multiplying F[0] by the encoded 
hologram, although the optical reimaging of the encoded symbols was 
not successful in this study. Unfortunately, recovery of the encoded 
information is extremely difficult due to critical alignment of the 
synthetic holograms. Recovery of the encoded symbols is possible, 
however, by use of the simulation routine and the results are plotted 
in Figure 40. 




Figure 32 . Simulated Autocorrelation of Input Character 
T with Matched Filter T. 



Figure 33* Optical Autocorrelation of Input Character 
T with Matched Filter T. 




Figure 3^. Simulated Output of Optical. Processor with 
Input Chaxacter and Filter Unmatched. 



Figure 35 - Optical Reconstruction of Autocorrelation 
with Unmatched Filter. 




Figure 36 . Simulated Autocorrelation of Input Character 
T with Aperture Position Shifted. 



Figure 37* Input Mask OXT for Character Recognition 
Study. 
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Figure 38 . Simulated Autocorrelation of Input OXT 
with Matched Filter T , 



Figure 39* Optical Autocorrelation of Input OXT 
with Matched Filter T • 




Figure 40. Simulated Reconstruction of the Hologram 
F[AB]/F[0]. 



CHAPTER V 


CONCLUSIONS AND RECOMMENDATIONS 
Conclusions 

The results of this research are oriented more toward appli- 
cations than to pure theory. Most of the conclusions drawn from the 
experiments are based on physical results. Of course, the mathematics 
of Fourier optics is the basis for the synthesis of Lohmann-type 
holograms. Also, techniques used for finding optimum signal to noise 
values are founded on mathematics. 

The synthetic holograms produced in this study reconstruct images 
as predicted by the theory of optical holography. The Gould plotted 
binary masks yield a hologram which produces a high quality image 
reconstruction. The optimization algorithm provides a method for 
determining parameters for synthesizing near-optimum holograms. 

Photographs of optical reconstructions reveal that the selected 
parameters produce holograms of quality as predicted. Figure 19 is an 
optical reconstruction of high quality as selected by the optimization 
routine. 

Selection of clipping and bandlimiting parameters may be made by 
varying those quantities and observing their effect on the resulting 
signal to noise ratio. Typical computer run time for the program SN 
to tabulate signal to noise ratios is in the order of 8 minutes. 
Application of non-linear programming techniques as in Subroutine OPT 
reduces computer time to approximately 3 minutes and directly provides 
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parameter values to the hologram plotting routine. Thus, program OPT 
saves both operator time and computer run time. 

Optical simulator plots reveal the effects of clipping and band- 
limiting variations on image reconstruction. The simulator- allows 
image reconstruction to be observed without the necessity of using a 
laser and optical bench. Note the high quality of reconstruction in 
Figure 23 with parameters set for optimum values. The simulated 
reconstruction of Figure 24 displays differentiated edges due to the 
clipping value selected. Finally, in Figure 25, the absence of high 
frequency terms due to bandlimiting causes a rounding of the leading 
edges of the image. For purposes of reconstructing images, the optical 
simulation routine provides reconstructions which agree precisely with 
optical reconstructions. 

Program RECON was devised to reconstruct the image from the 
hologram exactly as it was recorded by the plotter. The unique 
feature of EECON is that it takes into account both phase and ampli- 
tude quantization as introduced by the plotting equipment. The 
optimization program finds an optimum hologram in terms of clipping 
and bandlimiting only; thus, the important contribution of plotter 
induced errors has been disregarded in the making of the hologram. 

The plot shown in Figure 26 is a simulated reconstruction as recovered 
from the plot tape. The parameters for this hologram were the same as 
used in making the optimimi holograms of Figures 19 and 23. Obviously, 
there is more degradation than was present in the reconstruction 
through the simple optical simulator. 
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The additional degradation of the image reconstructed from the 
plotter is due primarily to phase and magnitude quantization as dis- 
cussed in Chapter II. In addition to quantization errors, there is 
some loss of magnitude and phase information due to occasional over- 
lapping of apertures in adjacent cells. A signal to noise study was 
run on the hologram recovered from the plotter. Main program PLOTSN 
w*as used. For the clipping and bandlimiting parameters selected, the 
signal to noise ratio of the hologram prior to plotting was 7.5454. 
After plotter degradation, the signal to noise ratio dropped to 3.3606, 
a decrease of image intensity to 48% of its original value. 

For testing the spatial filtering capability of the simulated 
optical processor, the Abbe-Porter mesh was applied as an input and 
a horizontal slit was used as a filter as shown in Figures 27 and 28, 
respectively. The result was a complete suppression of the horizontal 
components of the mesh as verified by Figure 29. Rotating the slit 
filter by 90® suppresses the vertical components of the mesh in the 
reconstruction as seen in Figure 30. 

One of the outstanding features of computer generated holograms 
is their use as spatial filters. Once a method for generating synthe- 
tic holograms has been devised, spatial filters for character recog- 
nition and autocorrelation studies may be synthesized with ease. The 
real advantage of synthetic binary spatial filter generation is that 
the necessity of using an optical bench, liquid gates, and recording 
apparati is alleviated. 

The program HOLPLO was used for generating synthetic binary 
spatial filters and the simulation routine SPAFIL was used for 
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simulating the results of the spatial filtering experiments. The 
experiments involved autocorrelation and character recognition 
studies. 

In the spatial filtering experiments, the symbol T was used as 
the filter. Various input characters were applied and the results were 
observed by both the simulator plots and the optical reconstructions. 
Figure 32 is a plot of the autocorrelation of the input symbol T with 
the matched filter T . Note that only the autocorrelation from the 
simulator is plotted. For comparative purposes, a synthetic binary 
filter of the symbol T was made utilizing HOLPLO. With a mask T as 
the input to an optical processor as sho^m in Figure 2, the optical 
reconstruction is photographed in Figure 33. In the photograph, the 
autocorrelation spot appears to the right, the input T is reimaged 
near the .center, and the convolution spot appears on the left side of 
photograph. Figures 34 and 35 show the results of an input character 
X which is not matched to the filter T . Note that there is some 
cross correlation, but not to the degree observed in the matched case. 
One further verification of the versatility of the optical simulator 
was shown by shifting the position of the s 3 nnbol T in the input aper- 
ture and observing the corresponding shift of the autocorrelation spot 
in Figure 36. 

A character recognition device was simulated by applying the input 
characters OXT to the optical processor. The simulated autocorre- 
lation with the filter T is plotted in Figure 38. Notice that the 
autocorrelation spot aligns with the character T of the input binary 
mask. A mask of the characters OXT was fabricated and used as the 
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optical input to the system of Figure 2 using the synthetic binary 
filter T. The optical results are photographed in Figure 39. The 
results obtained using the synthetic filter appear to be of quality 
comparable to those systems employing the Vander Lugt filters. 

One certain conclusion that can be drawn from the hologram 
coding study is that optical decoding of synthetic holograms is 
extremely difficult. Subroutine 0PT2 and a modified HOLPLO program 
were employed to find the optimum parameters and synthesize a binary 
hologram of the space domain division of two Fourier transfoirm arrays. 
Successful results were achieved employing the simulation routine as 
shown in Figure 40; however, optical recovery of the encoded symbol 
was never attained. The presence of the information in the encoded 
hologram was verified by autocorrelation, but due to optical alignment 
problems, a legible reconstruction was unsuccessful. 

There are two contributions of this study which are signficant. 
First, the optimization routine allows hologram parameters to be 
selected which produce the best possible computer generated hologram. 
The fact that the hologram is optimal has been verified by both 
optical simulation and by the laser and optical bench. The second 
contribution is the simulated optical processor which performed in the 
same manner as a classical optical processor. The simulation of image 
reconstruction allowed experimental results to be observed without the 
use of lasers and optical benches. The outstanding result obtained 
from the simulator was the image reconstruction directly from the 
plotter tape. Access to the plotted hologram allowed degradation 
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effects of phase and amplitude quantization to be observed as well 
as the standard computer constraints. 

Re commendations 

Based on the results of this research, it may be concluded that 
Lohmann-type holograms are feasible for use in data acquisition, 
storage, and transmission. There is definitely degradation of input 
information during the process of constructing the synthetic hologram, 
but in light of the relatively small amount of computer storage 
required to represent a large amount of data, some degradation is 
acceptable. Future studies should consider the use of Lee-type holo- 
grams and gray-scale plotters with an eye toward minimizing quanti- 
zation errors. 

A very promising area of application of computer-generated 
holograms is that of synthetic-aperture radar. Optical information 
processing techniques have been employed successfully in high- 
resolution synthetic aperture radar systems. Further investigations 
should include synthetic hologram recordings for these mappings. 

During the course of this study an attempt was made to observe 
contrast reversal effects in image reconstruction by filtering the 
zeroth order term of the Fourier transform hologram. This idea could 
be pursued with some possible applications in microscopy and image 
enhancement. 

The resxilts of the correlation and character recognition studies 
were as predicted by theory and verified by plots and photographs. 
Possible extensions of this phase of the study should include an 
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intensity measurement of the optical reconstructions and a computer 
algorithm for comparing relative intensities of the simulated- cor- 
relations . 

Finally, the computer routine for optimizing the synthetic 
holograms performed exceedingly well and produced optical and simulated 
results as expected. In the reconstruction routine RECON, there is a 
significant degradation of the reconstructed image due to plotter 
quantization effects. Therefore, it is recommended that a program be 
developed for determining optimum parameter values with plotter induced 
degradations being taken into account. 
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1* Literature Survey 

In an effort to appraise the state of development of computer 
generated holograms and spatial filters > an extensive survey of 
holographic related literature was undertaken. The dates of the 
survey extend from the beginning of computer holography in the mid- 
1960^s to the present time. The majority of this literature was 
available in the library of the Mississippi State University, A 
bibliography of computer generated holography and related topics has 
been compiled and appears in Appendix A. The sources of the articles 
listed and the inclusive search dates are as follows: 

Applied Optics , January 1967 through April 1976, 

2, Applied Physics Letters ^ January 1966 through May 1976, 

3, . Dissertation Abstracts International » b. The Sciences and 

Engineering, January 1967 through April 1976. 

Journal of Research and Development , January 1967 through 
March 1976, 

5 , Institute of Electrical and Electronics Engineers Proceedings , 
January 1967 through December 1975. 

6 , Institute of Electrical and Electronics Engineers Spectrum , 
January 1967 through May 1976. 

Institute of Electrical and Electronics Engineers Trans- 
actions on Computers , January 1967 through April 1976. 


7, 
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8. Journal of the Optical Society of America , January 1967 
through April 1976. 

9. Laser Abstracts , January 1972 through May 14, 1976. 

10. Laser Focus , January 1967 through May 1976. 

11. Nature , January 1967 through February 1976. 

12. Optica Acta , January 1967 through March 1976. 

13. Optical and Quantxim Electronics , January 1975 through May 
1976. 

14. Optical Engineering , January 1972 through February 1976. 

15. Optical Scanning News , January 1974 through May 1976. 

16. Optics and Spectroscopy , January 1967 through November 1975. 

17. Optics Communications , January 1969 through April 1976. 

18. Opto-Electronics , January 1969 through December 1974. 

19. Physics Letters , April 1967 through February 1976. 

20. Science, January 1967 through March 1976. 

21. SPIE Journal , January 1971 through December 1972. 

For the purpose of making a meaningful summary, the aforementioned 
topics will be broken into three categories : (1) computer generated 

holograms, (2) spatial filters and image processing, and (3) appli- 
cations of computer holography. Included in the summary of each 
category will be a list of the most pertinent articles with a 
snynosis . 


Computer Generated Holograms 

Several methods are presently in vogue for synthesizing holograms 


by means of a digital computer. Reference to a synthetic hologram 
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implies that a mathematical description has been used for forming an 
image and that a physical object need not be used. Of primary 
interest in this study is the binary hologram which has transmittance 
of either "1" or "0". 

The earliest references to computer-generated holograms date back 
to the mid-1960's -when Lohmann and Paris experimented with binary 
spatial filters. The binary Fraunhofer hologram synthesized by a 
digital computer was proposed in a paper by Lohmann and Paris. The 
hologram referred to in this article is known as the Lohmann-type 
hologram. This type of hologram subsequently has been improved in 
terms of reconstruction quality. Studies also- have been made of 
reconstruction errors due to quantization in Lohmann holograms. 

Two and three-dimensional holograms as well as color holograms 
have been synthesized by the digital computer. One of the recent 
synthetic holograms, is known as the Lee-type and employs the positions 
of the samples of the synthesized hologram- to record phase information 
of a complex wavefront* The kinoform is a computer generated wave- 
front reconstruction device which operates only on the phase of an 
incident wave.^^ More recently, a class of holograms called circular- 
carrier holograms (CCH) has appeared. 

Interest in computer generated holography appears to have 
revived recently as evidenced, by the number of articles and Ph.D. 
dissertations which have appeared during the last year. The current 
emphasis appears to be on applications of computer generated holo- 


grams. 
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Following is a chronological list of some of the most pertinent 
articles related to computer generated holography. A brief synopsis 
is included with each article in the list. 

1. Lohmann, A. W. and D. P. Paris, "Binary Fraunhofer Hologram, 
Generated by Computer," Applied Optics. Vol.6, No. 10, 1967. 

A method for snythesizing holograms from mathematical 
descriptions is described, A computer plots a drawing of 
the hologram which in turn is reduced photographically. 

Theory and experimental results are presented. 

2. Lesem, L. B., P. M. Hirsch, and J. A. Jordan, Jr., "The 

Kinoform: A New Wavefront Reconstruction Device," IBM 

Journal of Research and Development , Vol. 13, No. 2, 1969. 

A method for constructing an optical element which 
operates only on the phase of an incident wave is described. 
The kinoform exhibits high efficiency in terms of spatial 
frequency potential and reconstruction energy. Computer 
synthesis time is less than that for a digital hologram. 

3. Brown, B. and A. W. Lohmann, "Computer Generated Binary 
Holograms," IBM Journal of Research and Development . Vol. 13, 
No. 2, 1969. 

Improvements over the original Lohmann-type binary 
hologram are discussed. Some simplification in computer 
production and reconstruction of two and three-dimensional 
images are presented. 

4. Lee, W. H., "Sampled Fourier Transform Hologram Generated by 
Computer," Applied Optics , Vol. 9, No. 3, 1970. 

This paper describes a technique for determining a real 
non-negative function for representing the transmittance of 
a computer synthesized hologram. Positions of the samples in 
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the synthesized hologram record the phase information of 
a complex wavefront. This paper presents the Lee- type 
hologram. 

5. Gabel, R. A. and B. Lin, "Minimization of Reconstruction 
errors with Computer Generated Binary Holograms," Applied 
Optics , Vol. 9, No. 5, 1970. 

Degradation of image reconstruction in computer 
generated holograms due to equipment limitations is dis- 
cussed. Truncation and quantization errors are studied. 
Theoretical and experimental results are presented. 

6. Huang, T. S., "Digital Holography," Proceedings of the IEEE, 
Vol. 59, No. 9, 1971. 

This paper summarizes techniques for generating and 
reconstructing computer-generated holograms. Lohmann’s and 
Lee's methods of synthesizing holograms are included. Some 
applications of computer 'holography are presented. 

7. Chu, D. C. and J. R. Fienup, "Recent Approaches to Computer- 
Generated Holograms ," Optic^ Enjgineeri^, Vol. 13, No. 3, 
1974. 

A recent approach to computer synthesis of holograms and 
the production procedure is described. The method presented 
has the advantage of the kinoforra without its limitations. 
Experimental results are included and the referenceless on- 
axis complex hologram (ROACH) is introduced. 

8. Lee, W. H,, "Binary Synthetic Holograms," Applied Optics, 
Vol. 13, No. 7, 1974. 

This paper presents a method for making binary synthetic 
holograms of wavefronts with constant amplitude. 
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9. Lee, W. H. , "Circular Carrier Holograms," Journal of the 
Optical Society of America , Vol. 65, No. 5, 1975. 

Herein is presented a description of a computer generated 
hologram that is made with a reference wave having a linear 
phase variation in the radial direction. ' A method for making 
the circular- carrier hologram (CCH) and reconstructing the 
object wave front using a circular grating is formulated. 

10. Hugonin, J. P. and P. Chavel, "A Complement to the Theory of 
Lohmann-type Computer Holograms," Optics Communications , 

Vol. 16, No. 3, 1976. 

Brown and Lohmann improved their original holograms by 
using true phase and amplitude coding at the center of the 
diffracting aperture. This article presents a mathematical • 
analysis of this procedure. 

In addition to the ten articles listed above, the following 
references to the bibliography in Appendix A are also related to 


computer generated hologram; 

1, 

3, 

6, 

11, 

12, 

15, 

18, 

19, 

20, 

23, 

24, 

25, 

31, 

32, 

33-, 

34, 

35, 

36, 

37, 

38, 

39, 

41, 

42, 

45, 

48, 

49, 

50, 

51, 

52, 

55, 

56, 

58, 

60, 

62, 

64, 

65, 

66, 

68, 

69, 

70, 

71, 

72, 

74, 

75, 

76, 

77, 

79, 

81, 

82, 

83, 

86, 

87, 

89, 

90, 

93, 

94, 

95, 

96, 

98, 

99, 

100, 

101, 

102, 

103, 

105, 

106, 

107, 

109, 

no. 

111, 

113, 

116, 

117, 

118, 

119, 

120, 

121, 

122, 

124, 

127, 

130, 

131, 

133, 

134, 

135,. 

137, 

140, 

142, 

143, 

145, 

149, 

151, 

152, 

153, 

154, 

155, 

156, 

157, 

158, 

159. 








Spatial Filters and Optical Image Processing 

Computer-generated holograms are particularly well suited to 
spatial filtering applications because of their relative ease of 
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fabrication. In fact, filter functions which are difficult, if not 
dmpossihle, to produce by classical optical methods may be synthe- 
sized by digital computer methods. Furthermore, as is the case with 
computer -generated holograms, no physical object is required for 
synthesizing the filter. 

The development of computer-generated spatial filters closely 
parallels that of the computer-generated holograms. In 1966 Brown and 
Lohmann described a method for generating spatial filters by means of 
a digital computer. Their technique yielded a binary matched filter 
which exhibited the same, properties as an optically fabricated 
filter. 

In 1967 Burch Introduced a computer algorithm for synthesizing 
computer-generated holograms to be used as spatial filters.^ Lohmann 
and Paris proposed the use of computer created binary spatial filters 
for applications in coherent optical data processing systems. Included 
in their optical proc essin g applications were phase contrast demon- 
strations and gradient correlation filters. Another application of 
the binary spatial filter was its use as an inverse filter in an image 
restoration system. 

Several variations of Lohmann *s original filter have appeared in 
recent years. One of these is the kinoform applied to an incoherent 
optical processing system. The possibility of optical processing 
with incoherent light increases the scope of applications of computer 
generated spatial filters . Incoherent optical-image processing with 
synthetic holograms and the resulting signal- to-noise ratios has been 
studied in some detail. 
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r Computer-generated spatial filters may be synthesized in varying 
levels or shades of intensity as required for kinoform- generation. 

The use of a multiple gray-level plotting device for recording syn- 
thetic holograms has been reported by Campbell, Wecksung, and 
Mansfield. A halftone plotter has been used for image restoration 
in a computer simulated optical processing system. 

Following is a list of articles pertaining to computer-generated 
spatial filters and optical processing systems: 

1. Brown, B. and A. W. Lohmann, "Complex Spatial Filtering 
with Binary Masks," Applied Optics . Vol. 5, No. 6, 1966. 

This is one of the original papers on computer- 
generated holograms. A discussion of image restoration 
using binary holograms is included. Also of interest is a 
description of the use of a computer-generated hologram as 
a matched spatial filter. 

2. Burch, J. J., "A Computer Algorithm for the Synthesis of 
Spatial Frequency Filters," Proceedings of the IEEE, Vol. 55, 
No. 4, 1967. 

This letter describes a digital computer algorithm for 
synthesizing Fourier transform holograms to be used as 
spatial filters. Experimental results of filters made by 
this method are discussed. The Burch-type hologram is 
introduced . 

3. Lohmann, A. W. and D. P. Paris, "Computer-Generated Spatial 
Filters for Coherent Optical Data Processing," Applied 
Optics , Vol. 7, No. 4, 1968. 

In this paper the authors discuss binary spatial filters 
which are similar to Lohmann-type holograms. It is shown 
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that these filters can perform any optical processing 
operation producible by an optical filter; experimental 
results are included. 

4. Weider, H., R. V. Pole, and P. F. Heidrick, "Electron Beam 
Writing of Spatial Filters," IBM Journal of Research and 
Development s March 1969. 

A method for generating' spatial filters by use of 
electron beam--addressed KDP crystals is described in this 
article. ^ Some examples of simple filters made by this 
method are presented. 

5. ichioka, Y., M. Izumi, and T. Suzuki, "Halftone Plotter and 
its Applications to Digital Optical Information Processing," 
Applied Optics , Vol. 8, No. 12, 1969. 

A digital computer is used to simulate an optical pro-- 
cessing system. Examples of spatial filtering and image 
restoration utilizing computer holography are given. 

6. Patau, J. C., L. B. Lesem, P. M. Hirsch, and J. A. Jordan, 
"Incoherent Filtering Using Kinoforms," IBM Journal of 
Research and Development . Vol. 14, No. 5, 1970. 

Application of the kinoform as a filtering element is 
described in this article, and its advantage over the Loh-- 
mann type hologram is discussed. Some results of correlation 
filters are presented. 

7. Ransom, P. L. and R, M. Singleton, "Synthetic Real*“Imag inary 
Hologram Pair," Applied Optics , Vol. 13, No. 4, 1974. 

Two holograms, one representing the real part of the 
complex wavefront, the other representing the imaginary part, 
are synthesized. The computer generated hologram pair is 
illuminated and the resulting reconstructions are added in 
phase quadrature to produce the desired complex wavefront. 
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8. Campbell, K. , G. Wecksung, and C. Mansfield, "Spatial 
Filtering by Digital Holography," Optical Engineering , 

Vol. 13, No. 3, 1974. 

A multiple gray-level plotting device for constructing 
digital holographic spatial filters is presented. Examples 
of the synthesized filters in a coherent optical processing 
system are shown. The theory of sampled holograms and 
spatial filters is discussed in some detail. 

9. Lin, B. and N. C. Gallagher, "Optimum Fourier-Transform 
Division Filters with Magnitude Constraint," Journal of the 
Optical Society of America . Vol. 64, No. 9, 1974. 

This paper includes a discussion of Fourier transforms. 
The optimum spatial filter is developed using integrated 
squared error as a fidelity criterion with the magnitude of 
the filter transfer function subject to a constraint. Pri- 
mary considerations in this article are given to image 
transformation and various methods of improving image 
reconstruction. 

10. Chavel, P. and S. Lowenthal, "A Method of Incoherent Optical 
Image Processing Using Synthetic Holograms," Journal of the 
Optical Society of America , Vol. 66, No. 1, 1976. 

A Comparison of signal to noise in an incoherent optical 
processing system to that under coherent illumination is 
made. Image deblurring is discvissed, and image processing 
employing incoherent illumination on computer-generated 
holograms of a low number of cells is presented. 

Other papers, articles, letters, or books which are related to 
patial filtering and/or optical image processing are listed in the 
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bibliography of Appendix A. The following list of numbers indicate 


those 

entries in the bibliography: 

1, 

2, 

•4, 

7, 

9, 

10., 

14, 

16, 

20, 

22, 27, 43, 44, 45, 

46, 

50, 

53, 

59, 

61, 

73, 

75, 

81, 

84, 

98, 103, 104, 106, 107, 

no. 

111, 

115, 

119, 

129, 

131, 

132, 


133, 136, 146, 152, 154, 156, 165. 

Applications of Computer-Generated Holography 

The appeal of the computer-generated hologram is its construction 
from a mathematical description of any object — one, two, or three- 
dimensional. There are numerous applications of these easily fabri- 
cated devices. Among the list of applications are mass data storage, 
image manipulation, and optical element testing. 

An early application of holography was proposed by Gabor in the 
area of character' recognition.^ Later, this idea was extended to 
include Lohmann— type holograms and the recovery of coded information. ® 
Basically, these techniques are applications of the spatial filters 
discussed earlier in this paper. 

A promising utilization of computer-generated holograms is mass 
data storage of digital computers . Holographic storage offers the 
advantage of rapid accessibility and low cost compared with presently 
used mass storage media. 

Spatial filters generated by the computer have been used in 
optical image processing systems. Holography has been applied to 
the focusing of blurred images and to the enhancement of images 
degraded by instrument errors or other causes. This area of 
application is related to matched filtering and has been used with 
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both coherent and incoherent illumination. Other useful areas in 
computer holography include the reproduction of three-dimensional 
objects and color images. 

A very practical use of computer-generated holograms has been ■ 
the testing of optical surfaces. Several papers described aspheric 
lens testing and applications of computer holography to curved sur- 
faces . A related area of application is shearing 

interf erometry. In addition to optical testing, computer-holograms 
have been used for producing kinoforra lenses, zone plates, and 
diffraction gratings. 

Finally, a number of holographic apparati have been developed and 
patented. These devices include holographic recorders and character 
recognition machines. See U.S. Patent descriptions listed in 
Appendix A. The following list of articles applies to applications 
of computer-generated holograms: 

1. Gabor, D., "Character Recognition by Holography," Mature, 

Voi. 208, Oct. 30, 1965. 

A utilization of holography in the recognition of 
characters with many variants is presented in this article 
by Dennis Gabor, one of the pioneers in holography. It is 
proposed that holograms have properties which can discrimi- 
nate between ntnnerals and letters of the alphabet. 

2. Lohmann, A. W., D. P. Paris, and H. W. Werlich, "A Computer- 
Generated Spatial Filter, Applied to Code Translation," 
Applied Optics , Vol. 6, No. 6, 1967. 

An extension of Gabor's holographic character recognition 
systan is proposed which utilizes computer-generated 
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holograms. Some of the problems of producing the spatial 
filters are reported and experimental results are shown. 

3. Stroke, G. W. and R. G, Zech, "A Posteriori Image-Correcting 
’Deconvolution’ by Holographic Fourier-Transform Division,” 
Physics Letters , Vol. 25A, No. 2, 31 July 1967. 

Sharpening of defocused images and correction of instru- 
ment induced faults in photographic recordings are shown to 
be correctable by holographic methods. Although techniques 
formulated in this letter refer to classical holography they 
are adaptable to computer generated filters. 

4, Hickling, R., "Scattering of Light by Spherical Liquid 
Droplets Using Computer Synthesized Holograms,” Journal of 
the Optical Society of America , Vol. 58, No. 4, 1968. 

A proposal for determining diameters of spherical liquid 
droplets- in the range of 0.5 to 20 microns is presented. The 
method is based on techniques used in the construction of 
holograms synthesized from far field illumination scattered 
by the droplets. 

5. Jordan, J. A., P. M, Hirsch, L. B, Lesem, and D.L. Van Rooy, 
’’Kinoform Lenses,” Applied Optics , Vol. 9, No, 8, 1970. 

The production of accurate, thin, and light weight 
lenses is discussed in this paper. Large lenses which are 
very thin and light-weight compared with conventional lenses 
are constructed by kinoform techniques. 

6, Fillmore, G. Js., ’’Kinoform Viewed as a Random Number Gene- 
rator,” Applied Optics . Vol. 11, No. 10, 1972. 

Kinoforms are shown to produce output intensities that 


are a random selection from some universal population. 
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7. Engel, A, and G. Herzlnger, "Computer Drawn Modulated Zone 
Plates," Applied Optics , Vol. 12, No, 3, 1973. 

A method of producing modulated zone plates by means of 
a digital computer is described. The modulated zone plates 
are compared to synthetic binary holograms. 

8. Lee, W. H., "Contour Map Display Using Computer-Generated 
Holograms," Applied Optics, Vol. 14, No. 10, 1975. 

This paper presents an .application of Lee-type holograms 
in displaying contour maps of two-dimensional functions. 
Interesting maps of the force-lines of an electric dipole 
are shown. 

9. Takahaski, T., K. Kono, M. Kawai, and M. Isshiki, "Computer- 
Generated Holograms for Testing Aspheric Lenses," Applied 
Optics . Vol. 15, No. 2, 1976, 

A computer-generated hologram named the "inclined-bar" 
type is proposed for measuring aspheric surfaces. Examples 
are presented. 

10. Sirohi, R. S., H. Blume, and K. J. Rosenbruch, "Optical 
Testing Using Synthetic Holograms," Optica Acta, Vol. 23, 

No, 3, 1976. 

A modification of the Birch-Green computer generated 
hologram is presented. The ability of these holograms for 
evaluating aspherics is demonstrated., 

Refer to the following bibliographic entries of Appendix A for 
additional articles on the applications of computer— generated holograms 
8, 21, 29, 39, 40, 46, 47, 50, 54, 57, 63, 67, 71, 73, 


78, 85, 88, 91, 92, 108, 112, 114, 120, 123, 128, 129, 136, 138, 
141, 144, 147, 150, 162, 163, 164, 165, 167. 
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2 i A Compiled Bibliography on Computer- 
Generated Holography and Related Topics 


The following bibliography concerns the subjects of computer- 
generated holography, spatial filtering, and applications. -Some of 
the articles pertain to classical holography but appear to be appli- 
cable to computer-generated holography. 

The bibliography is divided into three sections as follows: 

(1) articles published in journals, (2) books, and (3) U, S. patents. 
The entries are listed chronologically except when the month and year 
of some publications are the same in which case the listing is 
alphabetical. An alphabetical author’s index follows the publications 
listing with references to the numbered articles of the bibliography. 

Some of thebibliographical entries are followed by a statement 
indicating the source of an abstract on that article. 


Bibliography 


Articles 
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"Optical Image Synthesis (Complex Amplitude Addition and Sub- 
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Vol. 18, No. 2, pp. 116-117, August 1965. 
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pp. 422-423, October 30, 1965. 

3. Waters, J. P., "Holographic Image Synthesis Utilizing Theoretical 
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Filtering by Means of Holograms," Ph.D. Dissertation, Xerox 
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Mich., 1966. 

Dissertation Abstracts International , B, The Sciences and 
Engineering, Vol. 28, 1967, p. 1505B. 
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APPENDIX B 


^ LISTING OF FORTRAN PROGRAMS 



A Listing of FORTRAN Programs 


The following programs and subroutines listed in alphabetical 
order are described in Chapter III, For a complete explanation of 
the use of these programs, refer to Chapter III. 





Ears *****-^ 


1000 


5 

■2 

1001 


subroutine bITS(A»ClIP»NQ»N12»B) 
INTEGER BTh600TiNPNT2> . 

COMPLEX A ( <34^64) 

DATA NPN(l) ,NPN(2 )/o777777777777»0/ 
INTEGER QT(2»39) 

COMMON/QUAM/QT 

NT=0TC1»NQ) 

NC=0T(2»N0) 

ACON=:FLOAT ( nC > /CLIP 
ANC0N=FL0AT (NO /6. 283185308 
Do 1000I~1>^JO 
WR1T£(9)(NPn<N12) »J=1»23) 

CONTINUE 
.NL=NC*23 
N2=NC/2 • 

CALL INIT(N2fNC) 

DO 2 I=lrNT 
DO 3 J=lfNL 

CONTINUE 
DO ^ 0=1 » NT 

R^-S^MINllREALtAdfO) ) »CLIP) 
IMAG=R*ACon+. 5 
IF(lMAG.EQ,0)GOTp4 
ANG=A1MAG(A(I» J) ) 

I ANG=ANC0N+ANG+ . 5 

!aNG=:MOD,(IANG»NC) 

i8L=(0~i)*NC . 

lVJL=IBL/36 

lSH=IbL-lWL*36+IANG 

IT=ISH/36 

IWL=IWL+IT 

ISH=ISH-1T*36 

call BS(B»IwL»1MA6»ISH) 

Continue 

NB=-22 

DO 5 II=1*NC 

NB=MB+23 

NFI=NB+22 

B'?JJ)=XOR?b'?OJ) »NPN(N12)) 

,WR iTE^of ( D ( J J ) # J J=NB » NFI ) 

continue 
continue 
DO 1001 1=1, 


n: 


WRlTEtO (NPN^Nl2)f J=l»23) 
CONTINUE 


(ETURN 

:nd 
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BITS2 ****** 


subroutine BlTS2<A»CLiPfNQ,Ni2) 


1000 


1001 


integer B(4bOO>{NPN(2) 
DA?A^NPNClfWpN(2)/0777777777777iO/ 


integer QT(2»39) 

COMMON/QUAN/QT 

NT=QT(1»NQ) 

WRlTE(9)lNPNtNl2) » J=:l»23) 
CONTINUE . 

NL=NC*23 

&L^INIT(N2»NCJ 
DO 2 1=1 » NT 
DO 3 J=1»NL - 
Bt J)=0 
Continue 
Do 4 J=lfNT 

R=AMIN1<REAl(A(I»J)) »CUIP) 

IMAG=R*AC0N+?5 

IFnMAG.EQ.O)60T04 

ANG=AlMA6(A(I'f J) ) ■ 

IANG=ANC0M*aN6+-«5 ■ 

IANG=MOD(IAnO»NC) 

IBL=(J-1)*NC 

lWL=IBL/36 

lSHsibL~IWL*36+IANG 

lT=ISH/36 

NB=-22 

NFl=N0+22 
Do 7 JJ=NB#NFI 
B(JJ)=X0R(8(JJ)»NPN{N12)) 
CONTINUE 

WRITE(9) <B( JJ) rJJ=Na»NFI) 

CONTINUE 

CONTINUE 

w8iTE? 9) {hiPNltNia) » J=lr23) 

CONTINUE 

RETURN - 

END 


^)RIGIIS^AO PAGE T.' 
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4.*; + + *+ HS *‘)^^^** 


5U) 

IfUT4f 


BS+ 


AGN 


S(0) 

SEED 

N 

w 


AXR’S 


LA 


AMA^U 

Al*l 

LA 

Ar^»GEED 

LA^U 

A/;»0 

OSA 

A2f 0»A1 

DS 

i\:> t W 

LA 


SA 

Al.N 

J 


OS 

ACirRS 

DS 

A^JrRS+2 

05 

A10»RS+4 

LA 

Al»+lfXU 

LA 

A2»0fXll 

AA 

A2»A1 

DL 

Ah » W 

LA 

AcS. !=3»Xll 

DSL 

AIj 1 0 f A3 

LA 


LAfU 

A9r0 

OSA 

AO » 0 > A3 

OR 

A?»AO 

LA 

A10»A9 

LA 

A5f-*^2»XH 

LA 

A'»fN 

AMA 

A4fA5 

SSL 

A4,i 

MSi^U 

A4.23 

AA 

A2» A4 

LA»U 

A4»0 

LXl 

A2,23 

AArU 

A4 >1 

OH 

Ao $ 0 r A2 

SA 

A7»0f A2 

OR 

A«»1>A2 

SA 

A9r l»A2 

OR 

Aln»2f A2 

SA 

AU»2»*A2 

TLL 

A4»A5 

J 

AGn 

DL 

DL 

AofRS. 
AOf RS+2 

DL 

A10»RS+4 

J 

b»Xll 

+0‘^00000000000 

RES 

1 

RES 

2 

RLS 

6 

END 
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CRUPT ****♦♦ 


SUBROUTINE cRUPT(b*C*CLIP»N9»N ) 
Complex A{b4»6‘i| rU(64f6*t> »c(04»64) 
InTEGlR QTT2»39) 

COMMON/QUAN/QT 

NT=QTtl»NQ) 

KcoN^FLOATtNO/CLIP 
ANC0N=FL0AT{NC)/6» 283185308 
NL=NT/2+2 , 

NK=(N+N-NT)/2+1 


1 

J=I»N 


DO 

lF((lt'LTlNL.0R.I.eT.NR).AND.<U.LT.NL.0R.J.GT.NR).0R.(NT»6C.M)) 


|^J)=C0. »0») 


>6^T04 

4 RsAMINl (REAL(B(I» J> ) 'CLIP) 
ANG=AIMAC(B(I 'U) ). 


ANG)»R*SIN(AN&)) 

CALL^FFT2{CrN» 

DO 3 I=j-*N 

BHs|bEAL ( C C I r J) ) <'*2+AlMAG ( C 1 1 1 J) ) t-*? 

AB= AB+RE AL { A ( I ' J ) ) *S0RT ( BMS ) 

CONTINUE 

iRR=tp*0SQ-2.*AB)»sC+ASQ 

ERR=ASQ/ERf! 

return 

end 


CRUPT 2 ****** 

r************************************************ ****’*'****■** ******"’■**** 
SURROUT I NE cHUPT2 (B'C'CLIP# NQ' N) , 

COMPLEX AXft4'64) 'DT64'64) »L<84f64) 

INTEGLR 0T(2'39) 
common/quan/ot 

NT=0T(1»NQ1 

NC=QT(2'NO) 

AC0N=FL0AT(NC)/CL1P , 

AnC0N=FL0AT ( NC ) /6 . 283185308 
NL=NT/2+2 
NR=(N+N-MT)/2+1 
DO 1 I=1'N 

iR(1i';/lHKl.or.i.6T.nr).and. (j.lt.nl»or.j.gt.nr> .or. (NT.GE.M)) 

GOTO I 

4 R=AM1NI(REAL(B(I'J))»CLIP) 

IR=R>iACgN+,b 

R=FL0AT<1R)/AC0N 

AMv'5=A1MAG{H(I'J)) 

I ANG- AN6* ANC0N+ ♦ 5 
ANG=FL0AT ( I aNG ) / ANCoN 
C(l»J)=CMPLXtR*COS(ANG) »R*SIN(AN6) ) 

1 CONTINUE 
RETURN 

end 
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****** ERR ****** 


^8a?K:rA!§8!H?:8fiH5af:t'?4Si6« 

Integer ot(2*39)- 

COMMON/QUAN/OT 

NT=OTUfNQ) 

•NC=QT(2r»Q) 

acon=float<nc»/clip ,, ^ 

ANC0N=FL0AT(NC)/6. 283185308 
NL=NT/2+2 


NR=(N+N-NT)/2+1 

88 i J=i!M 


IFC (I.Lt!nl,OR,I.GT.NR) .and. W.LT,NL.0R.J,GT.NR) .or. (NT.GF.N) ) 

>m?3.=io.,c) 

GOTOl 

4 R=AMIN1(REAL(B(I»0» )»CLIP> 

IR=R*ACON+.5 

R=rLOAT{IR)/ACON, 

ANG=A1MA6{8(I»J) > 

IANG=AN6*ANC0N+»5 ^ „ 

AN6=FL0AT( IaNG) /ancon 

ca» J)=CMPLXtR*COSC ANG) »R*SIN(AN6) ) 

1 Continue 
call FFT2(C»N) 

BSQ=0* 

AB=0, 

DO 3 1=1 »N 

BMS-REAL(C(I»U) )**2+AlMA6(C(Ir J) )**2 
BSQ=BSO+BmS 

AB=AR+REAL{AlI»J) )*SQRT(BMS) 

3 CONTINUE 

ERR=? SC+BSQ-2 . *AB ) *sC+ASO 

ERR=AiQ/ERR 

RETURN 

END 


****** CODE 

i:nLm**********************************************************‘^********** 

C MAIN PROGRAM CODE ^ 

Crt ****************************************'*^*****' i ******'***************** 
COMPLEX A(64f64) »Bl6A^»64) fEtfc't'GA) »D<64»64) »XX 
RtAL TE164) 

mo55‘tiL6, , 

READ (bflOOl) (TE<J) »J=lf64) 

Bo 1000 J=l,64^. 

B{IrJ)=CMPLxtTE<J) »0») 

1000 CONTINUE, ,, ^ 

CALL PIC(o.l.t39f64) 

Do 1500 in, 64 


1500 

1001 


2000 


100 


D(I.J)=CMPLX(TT<J)f 0.) 
CONTINUE, , ^ 

call PlC,(D{l.»39r64) 
F0RMAT<64F1.0), , , 
CALL ROTATE(Bf Af64) 
CALL ROTATE (0»E»64) 
CALL FFT2(A,64) 

CALL FFT2(E»64) 

DO 2000 1=1,64 
DO 2000 J=l,64 .. 
Xx=ATl»J)-E(IfJ) 
A(I»J)=XX 
CONTINUE 

call P0LAR(A»B»64)^ 
CALL R0TATE(BrAr64) 

DO 1 I=l»3 

RLAD(6flOO) NQfPCT 

F0RMATa2,F5»2) 

CLIP=PCT*FMaX(A»64> , 

call 13UILdIA»B»N0*64} 

call BIT5(B,CLIP*NQf2) 

CALL PIC(B»cLlPfNQ»6‘n 

CONTINUE 

EnU FILE 9 

STOP 

END 
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FMAX 

TUMCTION FmaX(A;N) 

COMPLLX A(6^»64) . 

SUM=0. . 

DO 1 1=1 »N 

I?(REAUAU,J)).GT,SUH)SUM=REAL(A(IfJ)) 
I CONTINUE 
FmAX=SUM 
RETURN 
END 


HOLPLO 

C***)k*3Si*^*****i*****«*=»!*.****************<‘*»'#*#****# 

c main program HOLPLO ^ ^ 

C***4'** «*♦****’*=■*=*= 4c ’»'**♦’''♦♦*♦**♦♦** 

COMPLEX A (64»64) 

-REAL TE164) 
oteNSION C(4600) 

DO .1000. 1=1,64 , , 

read (5t 1001 ) (TEtJ) »J=1»64) 


1001 


lOOO 


999 


loo 


98 

99 


g8'^SgJ^"3Ei:gl 

D(I» J)=CMPLX<TE(J) »0.) 
CONTINUE, , ^ ^ 

CALL PIC(D»i.^39f64) 
CALL R0TATEtB»Af64) 
WRITE(6»9a) , 

call FFT2(A,64) 

DO 999 1=1,64 

Do 999 0=1,54 
^{I,vJ)=CONJGUtI»J} ) 
CALL P0LAR(a»8,64) 

CALL R0TATE(B,A,64) 

Do 1 1=1,3 
READ{5,1Q0)NQ»PCT 
FORMAT C 12, F5. 2) 
CLIP=PCT*FMAX(A,64) 
CALL byiLn(A,B,NQ,64) 

loRpAlUUx, 2(311. 4>) 

END FILE 9 

STOP 

END 
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FFT2 ****** 


91 

20 


30 


100 

200 


1020 


1030 


1100 

1200 


SUF3ROUTINE FFT2(A»N) 
iNTEOtR gamma »5»6rU»D ^ . 

COMPLEX A(64»64) {bf Wt64WZ(64> 
GAMMA=AL0 G (float <N) ) /•6931471805+*! 

PHI=6»283185308/F'L0ATtN> 

bo 91 vi~l»M , , , 

THFT=PHI*FLoAT( J-1) , 

W( J)=CMPLX(&SCTHET> r-SlN(TH£T> ) 
CONTINUE^ , 

DO 200 IRslfN 
DO 20 1=1 fN 
7(1)=A(|»IR) 

CONTINUE 

S=1 

D=M 

DO 100 J=lr gamma 
DO 30 L1=1»S 
L1C0N=(L1-1)*D 
DO 30 L=1»D 
I=L+LiCON 

K=MOD*(*L+L1CoN*2 t N ) 

U=K+D 

iik}S?5IS}6, 

CorJTINUE 

S=S»2 

D=D/2 

DO 100 JJ=lrN 
Z<JJT=A<JJ,IR) 

CONTINUE 

CONTINUE 

M=N/2 

DO 1200 IR=1/.N 
DO 1025 I=1»N 
Z(I)=A(IR»I) 

CONTINUE 

Do'^liOO J=l, GAMMA 

DO 1030 L1=1»S 

I=L+UiCON 

IM=I+M 

K=M0D(L+LlC0N*2rN) 

U~K+D 

6=L1C0N+1 

b=2(U)*W(6J 

AaR»l)=Z(K:>t8„ 

AaRrIMJ=Z(K)-B 

CONTINUE 

S=S’3'2 

Bo^llOO JJ=1»N 

CONTINUE 

RETURN 

FnD 
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OPT •*=♦♦**♦ 


A^REFtSE FFT(FFT{A)) » WHILE 8 IS FFT(R) IN POI.AR FORM. 

COMPLEX A 1 64 » 64 ) » 8(64*64) 


real F(13*9) . 

CL=FMAX(B*64)/9. 

ASQ=0 • . 

DO 1 1=1*64 

ASO= ASq'+**rEAL(A(I* J) )**2 
CONTINUE 
N0=39 
l=NQ/3 
2 = 0 . 
x=o. 

tol=o.ooi 

KT=0 

30 CONTINUE 
Do 5 J=l*9 
NQ= 3*1 

fH »j)=err^aIb*asg*clip*nq,6'( J 

IF(F(I*J) ,LE*X) GO to 5 . 

x = Fy * J ) 

|K^BS(X-*Z) .6T.TOL) GO TO 5 

Continue 

WR1TE(6*110) NQrJl*X 
. J=J1- 

Do 10 1=1*13 
NCl=3*l 

CLIP=CL*FL0AT ( J) 

F(I* J)=ERR(A*B*ASQ*CLIP*NQ*64) 
IptFUfJ) .Le*Z) go to 10 
Z = Fd * U ) 

11=1 

ifIabsH-zi.gt.tol) go to io 

60 TO 25 


10 

CONTINUE 
WRITE (6*1: 
1=11 

25 


bO 

IF(<J1.EQ 
GO TO 40 

.55 

IFdl.EQ. 

11=11-1 

kt=kt+i 

^45 

GO TO bO 
11=11+1 • 
kt=kt+i , 

55 

GO, TO bO 
iF(Jl*tQ. 
J1=J1-1 
kt=ki+i„ 
GO TO '40 

60 

J1=J1+1 

kt=kt+i 

40 

N=im 

M=I1-1 


70 


65 

15 






KsJl-U , , 

Do 20 J=L*K 

IF((I.EQ.I1).AN0.(J.EQ.J1) ) go to 20 
F? iT5*)‘=ERR { a * B » ASQ * CLIP * NQ * 64 ) 

X=F(I*J) 

11=1 
J1=J 




|O^TO 25 


.. (KT.GE.3) GO TO 15 
20 CONTINUE 
CONTINUE 
WRlTCt6*99) 


Sl7i 



108 


OPT ♦♦♦*♦♦ 


99 


110 

100 

105 

115 


F0RMAT(///M the optimum bandwidth and clipping levels 

^WR1TE{6»100) li»Jl 

Nq=NQ-3_ 




^OrSaT t2xS^NQ= ’ » * C«-= » ► 12 { 3X » * S/N= ’ » G 12 . 6 , 

Format i6x» »nQ» *6Xf 'clip* »9x» !F ^ » I2» * » * » X2» * ) * » 

P8WSIiS;?f=^^{i?Sx:^3in!2nx.'S/N=',G12.6,//) 


//) 

//) 


RETURN 
END 


OPT! ♦♦♦♦** 


1001 

1000 


COMPLEX A (64 » 64) »B(64r64) 

gist 

RiAD?5jlOoi)^tTE(J> » J=l»64) 
F0RMAT<64F1.0) 

A?IrS)= CMPLX^TE(J) rO. ) 

•8^Ai!i!:^P^?2(A.64) , 

CALL polar (a»A»64) ^ , 

call opt(a*b»f»nq»pct»x) 

STOP 

END. 





0PT2 


c main program 0PT2 

COMPLCX A(64r64)fBU>tf#6^4) fCONJG 
COMPLUX E(6ii.64)»U«6‘*n>4)»XX 
RilAL TT(64) 
real TE{64) 

REAL F(13r9) 

DO lOOO I=l»64 

iooi mwusfiioS^^''”'''^^'*"' 

A?lrJ)= CMPLxfrEtv)) .0.) 
iOOO CONTINUE 


1500 


2000 


bo'Sk‘S?{ eh''" '''=""■" 
EiSfyApr'-’''”'-”-”-’ 
ptt mutiv, 

DO 2000 1=1,64 
Do 2000 J=i,64 
rr=real(e(i; j7) 

PP=“AIMAG(E{ I» J) ) 
b(I»jy=CMpLX‘RR»PP) 

CONTINUE 

call P0LAR{a*0»64) 
CL=FHAX(B,64)/9. 
call FFT2{a764) 

CALL F0LAR(a*A»64) 

ASQ=0» 

Do 1 I=l»64 
Do 1 J=lf64 

EEflfiffi * REAUA(1,J.|..2 
NQ=:39 
UNQ/3 
2 = 0 . 

X=0. ‘ 

TOL=0.001 

NQ= 3*1 

CLIP=CL*FL0AT{J> 

X=h { I » J) 

U1=J 

i^‘^8^>«x~z).gt.tol) go to S 

CONTINUE 

''^S1,T^**^'^10) NQ»J1»X 

U — J 1 

DO 10 1=1,13 
N0=3*I 

GO TO 10 

Z=P 1 1 » 

11=1 • 


ORIGINMi PAGE IS 
OP POOR QUAIOT. 


|g^^p||~Z).6T.T0L) GO TO 10 

“ fipSuo, „0I,.1,Z 
Go. TO 30 





0PT2 




45 


55 


60 

40 


70 


65 

15 

99 


11 = 11+1 

KT=KT+1 

•EQ.l) 60 TO 60 
Jl=Jl~l 

SSi- 

KT=KT+1 

N=ll+1 

H=I1-1, 

DO 20 I=H»H 
L=J1-1 

00^20 jJ=L»K 

lHKf^GT^0?^G0‘^T0 70 
IFid.EO.ll) .AN0.<J.EQ.J1J) go TO 20 
■ “ ‘ iUE 


C^NTH^y 


ERR(A»B»ASp»CLIP»NQ»64) 
^■'^GO TO 65 


110 

lOQ 

105 

115 


.. .. - -Le*X) 

X=F(I»J) 

11=1 

J1=J 

60 TO 25 

IF(KT.GE.3) GO TO 15 
20 Continue 

CONTINUE 

FORMATl///!' THE OPTIMUM BANDWIDTH AND CLIPPING LEVEL VALUES 

WRlTE(6rl00) I1»J1 
Nq=N0'"3 

ClIP= FLOAT ( *J1) /float ( II ) 

WRITE<6»105) N0»CLIP»X , ^ 

Format (2x» »no=’ » i2»3Xf *cl=» »i2f3X» *s/n=* »gi2.6»//> 


•Format 1 6x r » NO • »6X» »clip» »9x» ’f(* r ia» • » * » i2» • ) * »//) 
F0i<MAT<6X»l2>6'X»G8.4'3X*6l2.6) 

Format ( 2X » *i=» »i2»3x» »j=*»i2»3 
STOP 
END 


p 3X» *S/N=* fG12.6»//) 


oaGINMJ + f ^ 

OF POOR 



Ill 


3tc:4t3)E3(c4:i)( 


Pic 

SUBROUtlNF! pIC(Af Cl.lP»NO»N) 

COMPLEX At64»64) 

integer IG( 64 >-^^ 
integer 0T<2»39) 

COMMON/QUAN/QT 

NT=QT(1»NQ) 

NC-QTl2fNQ) 

WRITE (6» 101 )NQ»NT{NC»CLIP 
F0RMAT< ’1» »3ltO»GlO,b) 

acon=floatu1c>/clip 

DO 1 1=1 ♦ NT 
Do 2 J=l>h4T 

IG( J)=AMlNl(REALCA(lr J) ) » CLIP) *AC0N+*5 
CONTINUE 

VmiTE(6»l00) tIG( j) » J=lrNT) 

FORMAT (» »»&4I2) 

CONTINUE 

RETURN 

END 


POLAR 


SUBROUTINE pOLAR(A»R»N) . 
COMPLEX A(64»64) »B(64»64) 
DO 1 I=1»N 
DO 1 U=lrN 
BR=CABS(A(I,J) ) 
IF(BR.GT.O.)GOT02 
ANG^O* 

GOTO 3 , 


CONTINUE 

RETURN 

END 


BUILD 


SUBROUTINE BUlLO(ArB»NQ>N) 
COMPLEX A(64»64) »B(64»64) 
INTEGER 0T(2r39) 

common/quan/qt 

NT=QT(1»NQ) 

NC= 0 T< 2 eNQ) 

NL=N/2-(NT+l)/2+l 

NR=N/2+NT/2 

^0 1 I=NL»NR 
1N=IN+1 

DO 1 J=NL»NR 

BnN*f JN)=A(I» J> 

CONTINUE 

RETURN 

END ORK 



112 


bO 

c:00 


50 


102 

103 


c . MAIli f^fiOGtUM Ki:con 

c + * 4.* 4 •'-it* 4 4s4‘ + * + 'jc^+jc«-*l'4!+!i*4i4<SL»*v’5f**'* (,* + ♦++ + 

C HOST tis: RUil Sr.Olif-.NTlALLY aftkk holplo 

COfPII-f;R(r!-'1=AI5S) 

Pi\i<A‘.-t fi:i». hC = j2»t4T=64fHIO=39 
ulN.PfJblOit (,,(NT>lir) 
lilTlTOLR F tHc»MT) »L(f|T»NT) 

CO.-PUt-X A(bii»o!t) MJt64i64) »l;<6'^>6'0 
KUFGr.R IN (23) 
r>l=?j. 1415926 

1 CCfITTNUe 
1)0 60 1-1 r 40 
HbAlXy) IN 
com'iiiU'^ 
t)0 10 IG=l»ilT 
M=NC-#riC 
M-N-(iiC-l) 

JX=0 

Uo 50 JR=M»N 
JX=JX+1^ 

Ou 50 JC=1 ,hT 
F ( JX » 3G ) —0 
vJX-(» 

Uu !55 JL-M»N 
NX-JX + 1 
KtlAD(y) IN 
IC=fJ 

NO 15 0-1 1 
NA-WC/2 
lb=(NC/2-i) 
lU-NC/2 
IbX=0 . 

UjX-HC-1 
15-NC/2, 
a=FLOATUiA) 

,XU\=0 • 

C A— 0 • 

Z=x-1. 

JJ-36/NC 
NO 15 II-l/JJ 
IC=IC+1^ 

If( IC.GT.NT) ,00 TO l5 
F ( OX f 1 C ) CFLN ( 15 r 1 » I M U ) > . 

IF^r (JX»IC) .NH.O) GO To 103 

F^OX^XC)=FLut 

IF (F(JX»IC) .IC.O) GO TO 102 
IG=IS+7 

F(OX> lC)=FLN(ISi Ir IN(0) ) 

Xs-IS-4 
6p-TO_^103 

E?lGr^C^'=e(lGf IC)+F<0X» IC) 

15=15 me 
0K=R.OAT<E(lGiIC)) 

IFtnx.EO.O) on TO Pn 

IFUU^^Eo'^iV'rAND. (XCK.ET. .0001) ) GO TO 21 
II-(X.NE.CT) GO TO 20 
21 CONTINUE • 
fiL=0 
NJ!=0 
XCK=0. 

' r,0 TO 16 

ifick.eq.o.) go to 14 

Nl =Ml ^l 

If-‘(ML.F.().fJA) GO TO 7 

14 COMflNUE , _ 

if(I').fo.iox) go to 5 • 

Il)=t»-1 

/CK-=CK. 

GO TO 4 ^ 

3 IF UJL.LE.MC/4) go to Z 

GO TO 11 

2 NL=0 

XCK=0. 

7 ConriNUE , ,, 

B Crv=FLUtIN,l, IlJtjn „ „ ,, 

IFt (XCK.EQ. l) .AND. ICK.EO.O.) ) GO TO 12 
IFtCK.EO.O.) GO TO 3 
NR=NR+.l 


CpNTIR‘?!S 


_ ID.EO.IUXJ GO TO 
IU=I0+1 
XCK=CK 
GO, TO .0. 


IF 

IF 


R.LT.NC/'}»,,NR=9 
' .eq.na) • ANO.I 


MR.EQ.O)) go to 16 


OEffiWAIiP^ 
0? ?OOR OUAIIi^ 



'RECON + 

GO TO 13 

Ifj ■ ANL=r-LOAT(i'!L) 

THc.T=3*RI/.-:-ANL*2*P1/NC 
\jl) To 17 

11 AlJL=KLOAT(HL> 

TlltTrAfll.*2+i>l/HC-RI/2 
GO TO 17 

12 ANH=Fl.OAT(hR) , „ 

rHLT=Pl/2*-AMR'i‘2»-PI/MC 

GO TO 17 

13 AMf<=Fl.OAT(li|<) . ^ , 

THfc.T=*»50*PI/10n-AhR*2*Pl/NC 
I PTANR • ♦ 0 ■ ) TH£ T *0 • 

17 G*i^XO?lcT = yHET 

lU=NC/2*-(!C + ll 

l!)X=Il'.-NC/2 + l 
lL(X=IU<;HC/2-l 
20 CoilTKlUt 

15 COtJTinUE 

vi PoKMAVVlXf >1-'’ » lX*H3lS) 

10 conriNUE 

Do 61 1=1 » 'Ml 

oi Rl,ADty)_lN 

'♦3 KrITE(6»*-uV'^IC{ T » J) « J=l»MT) 

OO ^4 1=1 »Mr 

44 WKlTElb'93) lO ( I » 0) > 'J=i*NT ) 

93 FORMA ft IX » r„NGLC’ » 1X»32F3. U 

DO 42 1=1 » NT 
00 42 J=lfllT 

D(r»j)=c:!}n.x<t:(i''Ji ) 

format ( IX, 2SF4.1) 

Do 40 1=1, IfT 
Do 40 J=l»NT 
Rt-Ea»0)>C0S(Gtl»J)) 
AlM=FtI>0) i-5}ll(G(I,xl) ) 
ljTl»J)=CKPLXtKEfAIM> 

COMTIUUE 

FoRMATtlX,’ E’»IX»64I1) 
FORHATtlX,' E'»iXt2l4T 
call KOI ATE(B,A,64) 

00 600 1=1,64 
RE7%D (4) (E( I, J) ,J=1»64) 

DO GOi I-J»G4 
DO 6fji -0=1,64 
xx=a(i,j)=!<i:(I/J) 

A(I,0)=XX 
call FFf2(A,64) 

Call poLAfttA,«,64) 

Call ftoTATEtn,A,64> 

N0=39, 
pCT=.25 

CLIP=PCT!,FkaX(A,64) 
call uUIi.0(A,B,Nfl»64) 

CALL U1T52(b,CLIP,M0'2) 
WRlTCtojOa) 

90 F0RMA|l»in 

Call pic(b,clip,nq,64) 

500 CONTINUE 
ENOFILE 9 
STOP 
END 


42^ 

107 


40 

91 

100 


600 

601 


*^:>^^** ROTATE *■»■•=♦*'* 

SUUROUTIME rOTATE(A,U,N) 
Complex a 1 64 ' g4 ) , u ( 64 , &4 ) 
MA=N/2-2 
00 1 1=1, M 
DO 1 0=1, N 
lM=MOOtI+tlA,N) + l 
JN=MOU(J^NA,N)+l 
HUM,JN)=A(1,J) 

1 CONTINUE 
RETURN 
end 



DIMENSION U(64> 

DIMENSION V(64) 

Complex d{ 64*64) »E(64»b4) 

COMPLEX' C{b4»f>4) 

COMPLEX A{b4»64 ) »Bt64»64) >C0NJ6 
REAL TEC64) , 

REAO?b?lO01)^TE{J) * J=l*64) 

B a ► J) = cmplx^te{ J) » 0 . > 

COI4TINUE 

P0RMATU4F1.0),^ ^ 

call PlC(Bfi.f39»64) 

Call rotateTb»a»64) 

■WR1TE(6»98) 
call FFT2(A,64},, ^ 

CALL POLAR(ArB»64) 

HEAD{5»100)NQ»PCT 
call ROTATE «0»D»64) 

CALL build (DfE»N0i64) 

DO 41 I=lf64 

??J^=^{EAl1e?I*J) ) 

WRITE^6^l07) (U(J) r J5S1*64) 

107 FOHMAtVix»» ANGLE’r32F3.l) 

CaLL”cRUPtVb»C?CL1P»NQ,64) - 
CALL POLAR tc I B»64) ^ 

CALL K0TATElBrAf64) 

PCT=1.- 

*^lIp=PCT*FMAX « A 1 64 ) 
call BlTS2(A*CLIPfNQr2) 

100 rOKHATl I2rFS.2) 

WRiTc(6;?e) 

98 Format ( m m 

99 FORMAT C 4 ( 6X » 2C 1 1 .4 > ) , ^ 
call PIC(A.cLlPrNQr64) 

1 CONTINUE 

’ . 

END 

SfJ *:#»***- 

C main PROGRAM SN 

C********’f**’>f’f******■^•^r^^^’fm^^’^*■¥**<fF^.1r*^****■>|:*^*^*^f*:*■*.*t1Hr*■^ +*n<j(i**-*****-ii- 

Complex A (64»64) rB(64»64) 
real TE(64) 

REAL SN(9) 

COMMOr-i/OUAl VQT ( 2 f 39 1 
DO 1000 1=1,64 
READIS/1001)TtE( J) fj=l,64) 

1001 FORMAT(64F1.0) 

Dp 1000 U=l-,64 
aTi»J)=CHPLXITE(J) » o.) 

1000 CONTINUE 

call FFT2(A,64) • 

, call P0LAR(a»D»64) 

' CL=FMAX(D,64)/9. 

• CALL FFT2(A,64) 
call P0LAR(A»A»64) 

ASO=0. 

00 1 1=1»64 
Do 1 U=l>64 

ASO=ASO+REALfAtI,J))**2 

I CONTINUE 

DO 2 NQ=3»39f3 

WRITE{6'lOi)NOfQT(l»NQ) *QT(2fNQ) 


101 


100 

2 


FOliHATt///* 
Do 3 IFliig 


r3I5) 


SN^I)=ERRVA^B»ASQ»CLlP»NQrb4) 

CONTINUE 

V/RlTEt6fl00) (SN(I) »I=1,9) 
FOf?MAT(»0*,9012.6) 

CONTINUE 

STOP 

END 


OBlGrlNAL PAGE IS 
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sJ-iPLor 

+ ^c + 4 :^ + :4c^:4^i.* + * + V^4‘;ic**4;/Jc'ic V A + *: X A ♦«:**>?<* >}t ♦' + ‘5' + ^ ♦** <^ *>:* 'i' + + ’» J*f 4f** + 

c main PPOGiiAM SMHLOT 

c f«uGT ul: nuN sf:ouenti-ally aktg:r holplo 

C(M-’S^ILl'R(hLD=Al5S> 

PAi^A.'iui'rR fic=ia>MT=tf^^MIQ=33 
UI/^^nsIOM G(MT»MT)^ ^ 

IjTEGr.P f frJcM^T) rE(MjrNT) ^ 

COMPLEX A (64^6^U rb^6'^/64Wu(64?64) 
iNfgOLH, Iij(23) 

Pi= 3 ,l 41 L 92 b 


60 

^00 


'3 


lb 


II 


12 


13 


17 


?P 

lb 

bb 

92 
10 

61 

44 

93 


107 


40 


999 


60 


MR=0 

(nR.LO. 


0) ) GO TO 16 


. CO.-TINUE 
DO 60 X = W40 
READ (9) IM 
COf.TlNUE 
. OO in IQ^ir^iT 

C » iC*" 1 ) 

JX:=0 

Du bn JR=m^n 

^^R~^'lR+ 1 

continue 

If ( to.ta^iux) GO TO 9 
IO=tD+1 
:<Uf<.^CK 

GO To 0 
IP (NR.LT.nC/4) 

1P( (nl.co^na) 

Go TO 13 
ANL^PLOATUIl) 

ThL‘Tr3^PT/^-A^|L♦2^PI/MC 
Oo To 17 
ArjL=FLOAT(NL) 

TMET=Ai JL i^2^f>I/MC-Pl/2 
GO TO 17 
ANR^FLOAT (liR) 

TNLT=PX/2+AnR*2*PI/nC 
OO TO 17 
ANf<=FLOAT(NR) 

THCr=4bQ=i^PI/180-ANfU2->PI/NC 
IKTANR.tOtO.) THET-0. 

OO TO 17 
GilG/IO-TfU T 
lD-NC/2+MC4*ri 
Ib=NC/2-'t+NC + II 
1bx = Il>-WC/2+1 
lUX^IO+MC/a-l 

S8f;f?spi 

COt'/riNUE. 

Fok-M AT ( 1 X » t F ' f IX • 23 1 5 ) 

CoilTINUE 
DO 61 I=l»40 
REA0(9) IN 
Do 44 I=lrNT 

write ( 6 ; 9 J) (G< I » J) , J= 1 »NT ) 

f OKf :AtTix » ’ AUdLE • » lx » 32F3 1 1 ) 

Do 42 1=1 »NF 
DO 42 J=lrNT 

Ojl/JIyCFFLXtCtlfO) ) 

hOUMAT<lXf26F4.1) 

CO 40 t=i»nt 
DO 4n 0=1>NT 
f{E=E{l^J)i^CoG(G(If J>) 

AIM=EII>. 0 ) <-;;,INCGll»J) ) 

0(I»J>=C!-PUx{RE*A1H) 

COiJTTNUE 

UO 999^1=1 #64 

read (9T (0Tl»J)*J=lf64) 

ColiriNUE 
RE«0(9)FH 
fiLAD{9) FD 
/.::|-0/rkAXTG,64) 

WklTE(0»Aa) Z»FO*FD 

•FMAX0=MF10;4,E,X»»F:.‘AXfl=» 

call FFT2(0f64) 
call polar (D,0f64) 

CLip=TMAX<|t,64) . 

Call POLARCp»B»NT) 

CALL ROTAtE(B»B»NT) 

A5l— 0 • 

DO 1001 1=1, NT 
Do 1001 J=1,NT 


SlSoE QUAura 





S^JPLOT 


bO 


i02 

10^ 




in 


7 

t3 


looi 


600 


JX-JX+1 

Jx=0 

DO 55 «JL=M,M 
JX=vJX+l 

R£1A0(9) In 

1C=0 

Uy 15 j=l#23 
NA=MC/2 
lM=(NC/2-l) 

I0=NC/2 
lljX^O 
IUX=NC-1 
li-MC/2 
X-F LOAT(NA) 

XCK-0. 

Ck-O. 

2=X-1. 

JJ=30/NC 
DO 15 1I=1,JJ 
1C=TC-1 1 

l(-Uf;.5T.;!n GO TO 1.5 
r- ( JX> lO=FL(i(l5f 1» iiltO) ) 

I(-TF^(wX» IC) .NE.O) eo To 103 

F ( Jxf I c; ) =n.[j ( Ts » 1 » I ' i < 'J > ) 

IF£F( JX> IC) .Nt.O) GO TO 102 
1 I s+ / 

MJX>XC)::FLU(ISfl.lM(J) ) 

Is— IS**^ 
eo TO 103 
J5=fS+3 
CONTINUE 

tUG;IC}=E;{lG,IC)+F(JXrIC> 

lEj— IS'i'NC 

0K7FL0AT(£(iG,lcn 
IF (OK. HO. 0) GO TO 20 
CI=r LuAT ( n-jL) 

IF((C1 .E('.zy.ANO.(XCK,l 
IF(X.NE.CT) GO TO 20 
CONTINUE 

Nk=0 

XCK=0. 

CONTINUE 

CK=FLUaDrl,TN(U)) 

ifcck?eq!oI^66'''to*i^^ 

NL=NL+1 

If-(NL.EO.NA) GO TO 7 
LONTINUE 

IFdn.EO.Ifix) GO TO 5 

iu-in~l 

XCK=CK 

GO TO 4 

IK(NL.t-E.NC/4) 60 TO 2 

GO TO 11 

f<l.=0 

XCK-0. 

CONTINUE 

Ck=FLU(I0,1,IN(J)) 

IF ( (XCK.Lq.x) .and. ( cK.tQ.O. ) ) GO TO 1? 
IF(CK.EQ.O.) GO TO 3 
A5 (v=ASQ+PEA|.{D(I*U) ) 

Continue ' 

N0=30 , 

SlLE’g''^^" SN=.dGl2.6.//) 

STOP 

end 


.LT..0001)) 00 TO 21 
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***^** . GPAFIL 


c I'Aitf f'i’ugkam sPAf-ii: 


1000 


1500 

1001 


09 


<5000 


loo 


98 


AKi; nifni{Yr 'U£AL» At'm' 64 X'oH. 

■'fE , complex conjugate of the filter 
COMPL l-.X -A( o4j 64) »bl64;64) »tf64f 64) »l)(6(v»64) r^i-.r. 

HtlAL rL164) 

REAL TT(G4) 

00 1000 1=1,64 

i<t.An f8»l00J) <TE(J) » J=lr64) 

00 1000 J=l,64 

LiU» J)=C|.ipLx(TE<J) » 0 .) 

COMmiUE, 

call PIC(0,i.»39»64) 

Du 1500 1=1,64 

RLAlKbf lOOI) (-TT(J) » J=l»64) 

00 1500 J=l,64 
D(I,J)=CHPLxnT(J) , 0 .) 

CONTINUE 

CALL l‘lC(n,i.,39fb4> 

FORMA r<84Fl,0) 

CALL K0TArE(D,A,64) 

CALL i\OTATh. (0,E,b4) 

Call fft2(a,64) 

N0=33 • 

CLiP=.77*f-..'AX(A,64) 
call P0LAP.(/»,R,64) 

Call CHUPr^;(B,A»CLlp»Mo»64) 

Call fft2(f:,64) 
call- polar (E.MJ' 64) 

CLiP=.77*Ff-'AX(E,64) 

Call lPU!'T2 ( D , E , cl I P , mo » 64 ) 

Do 99 1=1,64 


UO 99 d=l,64 



=COtJOO<E(I» J) ) 

Do 3000 1=1,64 

UU 3000 J=t,64 

A(I,J)=A(I,j)*e(I,J) ■ 

coiiriNut 

Call eft? (a, 64) 

call P0LAR(a,8,64) 

CALL K0TATc;(U,A,64) 
PCT=1 • , 

Nu=39 

FOR.MAI (I2,J-6.2) 
CLIP=l'CTiFXAX(A,64) 

Gall uits?.(a,clip,no,2) 

WR1TE(6,98) 

Format (»i» J 

call pic (A, clip, NQ, 64) 

Continue 

End FILE 9 

bfuP 

END 


****** tape *♦-,1.** 


99 


SUUROUTINE TAPE(L,$) 

INIEGLR L(U 

Data Mon/o7760QOOogooo/ 

l?LA0(9,EMn=o9) (L il y , I = i,?3J 

Lr23)=AMD(Li23),M0D 

return 

HLTURN 2 
EMU 







